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SYNTHESIS OF HIGH-SPEED AUTOMATIC CONTROL SYSTEM STRUCTURES 


M. V. Meerov 


(Moscow) 


A method of synthesizing structures of automatic control systems allowing, 
in principle, an arbitrary accuracy of reproduction is obtained. A method is cited 
for taking into account the limitations resulting from the nonlinearity of certain 
elements’ characteristics, 


In previously published works [1-3] rules were established for constructing structures of automatic control 
systems which are stable under an unbounded increase in gain, In work [4] there was partially considered the 
question of selecting such a control system structure that, within its framework, it would be possible to obtain the 

necessary region in which the real frequency char- 
P - ‘ ‘ A acteristic of the closed system would be positive. In 
Pd fn net 
Din Bip) D,() Duo tm the works cited, stabilizing devices were employed, 
> I> the transfer functions of which had no free terms in 
the numerator. There was a large gain at zero fre- 
‘ quency (in the steady-state), In the present work we 
investigate the synthesis of structures within whose 
yy t= 7 frameworks may be obtained a virtually arbitrary re- 
gion in which the real frequency characteristic is 
positive, This means that systems may be obtained 
with a virtually arbitrary frequency pass band. Such 
systems have large gains, not only at zero frequency, 
but in an entire region of frequencies, Stabilizing devices for such systems have transfer functions of the type 
u /F(p). 


We explain the possibility of using, as stabilizing devices, elements with transfer functions of the type 
u/F(p), where p is a constant quantity. Examples of such devices are aperiodic, integrating, oscillatory and 
other of the simplest elements. It is clear that there is no difficulty in implementing such simple devices tech- 
nically. It is the more important to establish the static and dynamic properties which a control system might 
acquire if such links were introduced as stabilizing devices, 




















Fig. 1 


Fig. 1 presents the block schematic of an automatic control system. The system consists of N elements 
with transfer functions Kj/D;(p), of which n are shunted by the stabilizing devices in the form of local feedback 
connections with transfer functions y{/Fi(p). We assume that the gains of the n elements shunted by the stabili- 
zing devices may be varied within broad limits and, in particular, that they may be indefinitely increased. We now 
find the transfer function of the system under consideration. 


On the basis of Fig. 1 and of the rules for constructing the transfer function for complex systems, we may 
write 
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II D; (P) a K; 
“ K, - My jmn+i m (P) 
a ‘+Iaoll am 
(p) = "x, (1) 
1+ a ry KS 
= K, - Hy jenpi ; (P) 
i 
tow ao 


After simplification, the transfer function takes the form 


n N n 
Ws I] *; []*#s@ 





K(p)= i =1 j ~n+1 i:=1 m (2) 
n N n N n 
Hor%e+ex) T Fo+T es TT Aw 

t—3 j=n+1 i=1 jmn+i imi 


The characteristic equation of such a system will be 


n N n N n 
TM %e%+ex) [] 30+] TL x; [aw=o- (3) 


tm j=n+1 f=} j=n+1 i=1 


First of all, we determine what properties the system will have in the steady state. For this, we substitute 
in (2) p = 0: 


K]] 
K (0) =— t=! 7 - , (4) 
Wi%O%0+%,%) T] 42O+KT] 4. oO 
foi 


j=n+1 t=-] 





n N 
where K= Il K, I] K, 1s the over-all gain of the open-loop system. 
—a wet 


For static systems I] D,(0) = 1. 


im] 
n 


Assuming that |] F, (0) = 1, we obtain 
imt 





K (0) = x 


5 
[] 1+ Key+ ” 


t-1 


It is clear from (5) that the larger 4, the larger will be the steady-state error. However, even in the 


given case, by choosing yj < 1 (and, if necessary, uj << 1), we may obtain the necessary degree of steady-state 
accuracy, 


We assume that Kj increases without bound. With this, the system will remain stable if the following con- 
ditions hold (1-3), 


Denoting the degree of Dj(p) by dj and the degree of Fi(p) by vj, we have, on the basis of previously 
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obtained relationships [1-3], dj = 2+ 0— yj since, in the given case, the power of p in the numerator of the 


stabilizing device's transfer function equals zero. It follows from this, as a structural property of systems of the 
type being investigated, that a necessary condition for stability as Kj -> o is 


dj + vi = 2. (6) 


Three cases are possible here; 
a) vj = 0. This is the case for a rigid feedback connection; with this, dj = 2. 


b) vy =1. The stabilizing device is an integrating or aperiodic link; with this, the shunted element may 
be of first order, since dj = 1 and vj = 1. 


c) vi = 2. The stabilizing link is second order, and only a noninertial amplifying link (dj = 0 and vj = 2) 
may be shunted. 


We investigate all three cases in more detail. 


a) Introduction of a rigid feedback connection, In this case the polynomial Dj(p) may be of zero, first 
or second degree. Substituting in (2) F{(p) = 1, we obtain the transfer function for the case under consideration: 








n N 
Ix I] *; 
Ka (p) - n <1 N yoni n N . (7) 
Teo+e) TT 2+ I %; 
i—1 j=n+1 im] jeen+1 
The characteristic equation is written in the form 
n N n N 
IT i+ Xie T] Dip+ TX TT Xi =0. (8) 
i—1 j=n+1 i—1 3 jen+i 


The D-subdivision curve for the over-all gain is defined by the equation 


n N 
K = — [|] (Di(jo)+ Kim [] Di(ja). (9) 
i=} jon+1 


Transfer function (7) falls under the symmetric case [4], and the nature of the system is completely deter- 
mined by the properties of the D-subdivision curve. 


In the given case, the degenerate characteristic equation, as Kj —* 00, will be 


n N N 
Tl v: I] Di(~)+ T] xs =0, (10) 
i—<1  jen-+i jon+i 


from which it is clear that system stability, as Kj > oo, will essentially depend on the magnitude of the pj and 
on the parameters of that portion of the contour which is not shunted by stabilizing devices. Moreover, it is 
necessary to require that the auxiliary equation of the first, second or third type [1-3, 5] (which depends on the 
degree of Dj(p)Jalso satisfies the stability conditions. 


We assume now that the parameters are so chosen that the system remains stable as Ki-> @, The nature 
of the system is completely determined by the degenerate transfer function 


613 











N 
Il ¥; 


j=n+1 





K deg (p) = n N 


He TL tp) + I K; 


i=. j=andi jon+i 


: (11) 


N 


obtained from (7) as Kj > o or, for selected Il K from the placement of the D-subdivision curve for 
j=n+1 


N n 
Kaeg = — J] Dio) J] ws. 


j=n+1 i—1 
Thus, it is possible to draw the following conclusions as to the properties of the resulting system: 


1. The steady-state accuracy depends on the over-all gain and on the coefficients u;, where the steady- 
state accuracy will be the higher, the larger Ky and Kj and the lower pj. 


2. The domain of positiveness is limited by the critical frequency for D-subdivision for the gain Kdeg: 
and it is not in any way possible to enlarge this domain of positiveness within the limits of the given structure. 
n N 
3. The magnitude of the overshoot is determined by the character of the curve Kdeg* I] ey I] D j(jw) 
i—1 9 j=n+1 
and the selected (within the stable region) magnitude of the degenerate system's gain. 


However, even in this simplest of cases, the dynamic properties of the system are essentially improved as 
compared to the original one-loop system without the rigid feedback connections. 


b. Stabilizing devices with transfer function of the type wi/(1+tjp). In this case, in order that the system 
remain stable as the gain increases without limit, the degree of Dj(p) must not be higher than one. 





We assume, for simplicity, that Dj(p) is also an aperiodic link with transfer function K;j/(1 + Tj p). 
According to (2), the transfer function for the given case will be 





Kp(p)= 5A (a2) 
where 
n N n 
A= f] xk, J] ¥;[[ @+%>). 
im] j=n+1 i=} 
n N 
B= [] (4+ 204+77)+Ki4) [] +7, 0, 
i=] j=n+1 
and the characteristic equation has the form 
n N 
Hit+apaQ+7ip)+Kivd T] 4470p) + 
imi jon+1 
(13) 


n N n 
+1] Xi I] xi[]G@+7)=0. 


imi j=nt1 im} 


N 


The D-subdivision curve for the over-all gain K = || K; is defined by the equation 
ad | 
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n N 
[] +7 feit+%70)+ Ku) T] @+7;7o) 
K = —.— —2z4 ; (14) 


[] +474) 


i—j 





n 


Dividing numerator and denominator of Expression (12) by II (1 + % p) and substituting jw for p, we get 
i—1 





. K 
K, (7 @) = ; N (15) 
[I] 47) 44+7,i04+4,4) JT] 44+7;7) 
K + i—] j=n+1 





n 
]] @+ 7) 
i=1 
It is clear from this that the dynamic and static properties of the system are completely determined by the 
character of the D-subdivision curve for the over-all gain (here we have to do with the case of a symmetric trans- 
fer function [4]). 


The question of stability for unbounded increase of the gains Ky will not be considered here, since the in- 
vestigation of stability may be easily carried out in the same manner as in [1-5]. 


We mention only that the basic stability condition, dj + vj = 2, holds in the given case (dj = 1 and wj = 1). 


For stability as Kj —> o, it is necessary so to choose the parameters of the stabilizing devices that the de- 
generate equation and the auxiliary equation of the second type individually satisfy the stability conditions, For 
Kj > o, the dynamic properties of the system will be completely determined by the degenerate transfer func- 
tion, which is ~btained from (12) by allowing Kj —> © in it, i.e., 


I chet wn 


j=n+1 i—1 





Ky deg(P) = = N (16) 
ITu TJ @+7;~+ Il KT] tts) 
i—1  j=n-+1 j=nti. = i= 


The dynamic properties of the degenerate system are completely defined by the character of the D-sub- 


n 
division curve for Kdeg = I] Ki, Indeed, dividing both sides of (16) by II (1 +- % p) and setting p = jw, 


j=n+ti1 i—1 
we get 


, K, 17 
Ky deg (J) = “Seah ad : pe? 


Ile I] ¢+7;7% 


K gegt t= _innt 
[] 4+ 7%) 


imi 








The denominator of (17) consists of the sum of Kgeg and the right member of the equation of the D-sub- 
division curve for Kdeg, and the numerator consists of ry quantity Kdeg (the symmetric case [4]), It is now easy 
to determine the Seat structural properties of the system. Besides the necessary steady-state accuracy (which 
remains a finite quantity), which can be obtained by the proper choice of wi, Kj and Kj, it is possible to obtain 
an arbitrary region of positiveness and the necessary limitation on overshoot, With this, the region of positive- 
ness, in the given case, is not limited,in principle,to the critical frequency wp of the D-subdivision curve for 
the gain for the elements not shunted by stabilizers. In fact, the equation of the D-subdivision curve for the gain 
of the elements not shunted by stabilizing devices is written in the form 
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N 
K = = I] 447570) =D, (jo). (18) 


The equation of the D-subdivision curve for the gain of the degenerate equation will be 





n N n 
I] I] (i+ 7; j) [] #: 2% Go) 
Reg hil ai i—] _i=att om sea , (19) 
I] @+ 72) I] @+%7 4) 


i—1 i—1 


It is clear from (19) that the introduction of aperiodic links as stabilizing devices, followed by an increase 
in the gains Kj, increases the region of positiveness, If the difference in the powers of w in the numerator and 
denominator of (19) equals two, then the region of positiveness may, in principle, be increased to infinity, since 
the difference in phase here, as w varies from zero to infinity, cannot exceed # . The number of aperiodic links 
which must be introduced into the system in order to guarantee such a phase difference is defined by the relation- 
ships 


N- 2n =2 or n= -1, (20) 


where N is the total number of aperiodic links in the original one-loop circuit. (if N is odd, thenn = 





-iJ 


Thus, the use of aperiodic links as stabilizing devices, with the conditions cited above, gives the possibility of 
increasing indefinitely the region of positiveness of the real frequency characteristic. 


N+1 
2 


It is also easily seen that the cited type of stabilizing device lowers the magnitude of the overshoot. In- 
deed, as follows from (19), with the introduction of the stabilizing devices, the curve Dy (jw) is rotated clock- 
wise and is reduced in scale. 


If the variation in scale is compensated for by the corresponding choice of Kgeg then, precisely by this, 
there will also be attained the necessary region of positiveness and the necessary limitation on the overshoot. 


It should be kept in mind that here, in contradistinction to the case when flexible feedback connections 
are introduced, the magnitude of the time constants rj of the stabilizing devices has a greater effect on the 
character of the control process. We shall not stop to consider the other combinations which are possible here, 
for example, the case when integrating links are used as stabilizing elements, or when a single-loop circuit con- 
tains n integrating links which are shunted by stabilizing devices with transfer functions of the type 


oo 
1+ ™|P 
In these cases, the qualitative results are similar to those obtained above, and in each individual case, 
using the method presented here, it is not difficult to obtain quantitative relationships as well. 
c. Stabilizing devices with transfer functions of the type ap +" “ye In this case, in order that the 


system remain stable under unbounded increase of the gains, the stabilizing devices of the given type must shunt 
only a noninertial amplifying link. 


We consider a control system which consists of N aperiodic links and contains n amplifiers, the gains of 
which can increase indefinitely. We shunt each of these amplifiers with a stabilizing device with the transfer 
function 





Bi 
ajp’ + bjp + 1 


The transfer function of such a system, in the closed state, may be obtained if in (2) we substitute (1+ Tip) 
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for Dy (p) and (ajp* + bjp + 1) for Fi(p). As a result of these substitutions, we obtain 


n N n 
II K,]] K,|]| (a, p? + & p + 1) 
K (P= = ti — (21) 





[] G@irtt+opt+i+ Ke) []G0+7;, 24+ ]] %[] XJ] Grt+oye+t 


im] j=1 i—1 j-1 t=] 


The degenerate transfer function for Kj -> o is written in the form 


N 
va | X; 
Kieg= n N ah . (22) 


T]eJ]@+%7) wy 
iw] j~1 +]] K, 


n 


Tl@et+he+s *™ 


i—j 








It is obvious from Expression (22) that the system's properties are completely determined by the character 
N 


of the D-subdivision curve for [|x ie 
j=1 


We now consider the characteristic equation obtained from the degenerate system 


n N N n 
I] wT] + 57) + 11%; T] ap? +o:p +1) =0. (23) 
1 -1 j=-1 i=—] 


Keeping in mind that the system is stable for Kj > oo [aj and by are so chosen that the auxiliary equation 
of the second type and the degenerate Eq. (23) individually satisfy the stability requirements], one may conclude 
that the derived degenerate system is equivalent to the original system at the input of which are applied 2n ideal 
derivatives, Here, the closeness of the derivatives to the ideal will be the greater, the larger the amplifier gains, 
Ki. In such a system it is possible to obtain an infinite region of positiveness. 


For this purpose the number of amplifiers and the corresponding number of stabilizing devices may be 
determined from the relationship N-— 2n = 2, or 


oe 0 


If the original system consists of N links with transfer functions Kj/ Dj(p), then it is proven, by considera- 
tions analogous to those given earlier, that an infinite region of positiveness may be guaranteed if the number 


of noninertial amplifiers, and the corresponding number of stabilizing devices, is defined by the following re- 
lationship: 


N 
Dd —2n = 2, (25) 
i—1 
where dj is the degree of the polynomial Dj (p). 
‘ N 
From this, we get that n = zm d,;—1. 
im 


The degree of overshoot may be limited by the proper choice of the parameters aj and bj, which follows 
N 
directly from the expression for the D-subdivision curve for |] K;, since 
j=l 
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N 
W ]] @4+ 7,7) 
K =||K; =— ~ ix 


St [] fa Go)? +570 +4) 


im 





(26) 


The method considered here for improving the dynamic properties of a control system can be realized by 
the artificial introduction into the system of an amplifier with a large gain, shunting the stabilizing devices of 
the type stated, or by the use of high-gain amplifiers which are already in the system. 


The structures obtained, as a matter of fact, render possible the obtaining of high gains, not just at zero 
frequency, which was the case with the introduction of flexible stabilizing devices into the system, but over the 
entire frequency range, to within sufficiently wide limits. 


Some General Considerations 





In order to present a more complete physical picture of the method given in this work of improving the 
dynamic properties of automatic control systems, we consider individually the portions of the circuit with high 
gains which are shunted by stabilizing devices with transfer functions of the type 




















a 
a i/ Fi(p). 

Zin > £ Fig. 2 gives the schematic of such a portion. In the direct circuit there is an 
element with the transfer function Kj/Dj(p) which is shunted by the stabilizing de- 
vice w4/F{(p). The transfer function for the circuit of Fig. 2 is written in the form 

Pp out (P) vad K;F ; (Pp) (27) 
: %n (P) Dj (Pp) Fy (Pp) +K yu, * 
Fig. 2 
Dividing numerator and denominator of (27) by Kj, and letting 1/ Kj = mj,we obtain 
our (P) _ F; (P) ; (28) 
Tn (Pp) m, D,(p) F {p)+ #4 


We now assume that Kj —> oo and, consequently, that mj—> 0. Then (28) degenerates to 


Zout(P) _ Fe(e) 
Zin(p) B; 





(29) 


If there are n such portions connected in series, and the gains of all these portions are sufficiently large 
(theoretically unbounded quantities), their over-all transfer function, for Kj» oo, degenerates to the following: 


n 
® out (P) ra F; (Pp) : (30) 
Tin(P) ing Hi 


It is clear from (29) and (30) that for unbounded increase of the gains Kj there is, on the one hand, an 
inversion of the transfer function  {/F{(p); depending on the character of Fj(p),there are introduced into the 
system ideal derivatives with their highest orders equal to the degrees of the polynomials F{(p) and, on the other 
hand, the influence of Dj(p) on the system dynamics is removed, However, such a system may be realized only 
if it remains stable for Ky oo. 


There is nothi:g new in such an explanation of the results obtained, In work [6] there is given an inter- 
pretation of the results of [1], starting from the ability of an amplifier with a feedback loop to invert the transfer 
function of this feedback path if the amplifier gain is sufficiently high. In work [7] it is also stated that the 
shunting of a high-gain portion of a circuit by a feedback path with transfer function P(p)/ Q(p) is equivalent to 
the introduction, in series with the original circuit, of an element with transfer function Q(p)/P(p). It is stated 








with this that, if the time constants of the stabilizing devices are chosen equal to the time constants of the 
original circuit elements, then the operators with the identical time constants are cancelled out. 


It is essential to emphasize that the choice of the time constants of the stabilizing devices is not arbitrary, 
even in those cases when the amplifier gains are as high as desired. The fact of the matter is that the stability 
conditions, either for the degenerate equation or for the auxiliary equation, essentially depend on the time con- 
stants of the stabilizing devices. If it is considered that the degenerate equation, by its dynamic properties, is 
equivalent to a system described by a first- or second-order equation, which is always the case for a system with 


an unbounded region of positiveness, it is then seen 











MK Ks Ng a; i Lt, that the realizability of the system depends com- 
I*hp M+], {+h,p = fe /+], /+i,p 





pletely on the stability condition for the auxiliary 
equation being satisfied. As an illustration, we give 
the following example. Fig. 3 gives the block sche- 











matic of an automatic control system, The initial 
tp te itp single-loop circuit consists of six aperiodic links. 
Three stabilizing devices, with transfer functions 
Fig. 3 ni /(1 + Tp ), are introduced into the system in 


order to obtain the required dynamic properties. 
In this case, the transfer function of the system is written in the following form: 


A 3 
K (p) = a4 (31) 


where > . 
A = K,K,K,K.K,K, || (1 + ™P), 


t—1 


8 
Ba=[](4+7:p)(1 + Tap) (A + up) + mK, 


im} 
C = ((1 + Tp) (1 + top) + uakKs] ((1 + Top) (1 + top) + peXs), 
3 
D = K,K; || (4 + up). 
i=—1 


We assume that the gains Ky, Kg and Kg are sufficiently large and, for simplicity, are all equal; K, = 


= Kg = Kg = K. We divide (31) by K® and let K—> oo; then the degenerate transfer equation takes the following 
form: 


3 
Ki KiKs II (i ao ™P) 





K,., (p) = ‘<1 ; 
de 3 3 
‘ vitatell (1+ Typ)+KikoKs I (1 + tp) (32) 


The degenerate characteristic equation will be 


i=—1 


3 3 
wits TL (1+ Tip) + KiK Ks [] (1 +p) =0. (33) 
t—1 


It is clear from (33) that there is no difficulty in stabilizing the degenerate system, In particular, if 
Ti = 74, then Eq. (33) is written in the form 


3 
Il (4 + Tip) (vitots + K,K,K5] = 0 


im] 


and the stability conditions are satisfied for any Tj, Kj and pj. The auxiliary equation is decisive, since on the 
satisfaction of the stability conditions for it depends the realizability of the system, We now set up the auxiliary 
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equation, We assume, for simplicity, that #, = fz = 3 = 1. After some simple computation, we obtain the 
characteristic equation of the complete system in the following form: 


6 3 3 

mI] (1 + Tip) II (1+ tip) +m*]] (1+ Tip) ((1 + Top) (1 + up) (1 +7 5p) x 

=1 =] =} 

x (1 + tp) +(1 + Tap) (1 + ap)(t1+Top)(1 + top) (1 +7 5p) (1 + tp) A+ 
3 (34) 

+T7,p)(1+sp))-+-m |] (4 + Tip) (4 + Tap) (1 + ap) + (1 + Typ) (1 + ep) + 


im} 


3 3 
+ (1 + Top) (1 +top)] +] (1 + Tip) + KK 2K [] (1 + cup) = 0(m=t/R). 


In the given case, as is obvious from (34), the auxiliary equation is of the second type. In order to find this 
equation, it is necessary to find the two leading terms in each of the polynomials multiplied, respectively, by 
m*, m*, m and m’ . We consider the case when all the time constants T andr are identical. For this case, 


the auxiliary equation takes the form 
Tg? + 67%q* + 37%q* + 127393 + 37%q? + 67g + 1+ K,K,K; = 0. (35) 


By setting up the Hurwitz determinant, we easily convince ourselves that the stability conditions are not 
satisfied for any value of T. Consequently, such a system cannot be realized. 


Thus, it is impossible to draw conclusions as to a system's properties based only on the fact that operator 
factors may be cancelled out due to the inversion of the transfer functions of the stabilizing devices shunting 
high-gain amplifiers, This can lead to serious errors. Therefore, in all cases, it is first necessary to ascertain 
the stability conditions on the auxiliary and degenerate equations, and only then, on the basis of these conditions, 
to choose the time constants of the stabilizing devices, It should be kept in mind that it is not necessary to 
strive for absolute accuracy in reproduction by equating the time constants of the stabilizing devices and the 
time constants of the links not shunted by stabilizing devices although, in the general case, if the time constants 
in the original circuit are not identical and py; # wy # ws # 1, such complete compensation is also possible. 


On Allowing for the Power Limitations of the Individual Links 





The practical realization of the structures obtained does not present difficulties, However, when this is 
done, there may arise inadmissibly large values of the variables (for example, currents in motors, etc.) in the 
individual elements, Some of these variables may enter regions where the characteristics of the elements be- 
come nonlinear (for example, engender saturation). The consideration of saturation amounts to this, that the 
averaged values of the amplifier gains Kj will have finite values, and mj = 1/Kj will also have finite values. 
This means that the derivatives will not be ideal and, naturally, that limits the upper bound on the frequency 
pass band. However, for practical usage, it is always possible so to choose the amplifier parameters that the 
averaged gains will not be less than some previously given values, so that the corresponding values mj = I/Kj lim 
will be less than some given value m,. Computation of the quantity m, can be carried out, for example, as is 
done for the quantitative determination of small parameters [5]. Significantly more complex are the questions 
related to the inordinate increase of certain coordinates, In order to remove such a possibility, it is necessary 
to introduce corresponding limiters into the system, It is for this purpose, for example, that cut-offs on current, 
voltage, etc.,are introduced into electro-drive designs. As a rule, such limitation occurs only in certain elements. 


In the theory of optimal systems [8,9], these elements are separated out as the so-called "nonvarying” 
portion of the system. For this portion of the system, which is described by an equation of lower order than for 
the entire system, the entire arsenal of the theory of optimal control may be employed. If use is made of a 
method for determining the characteristics of a nonlinear link designed for coordinate limitation, as this is done 
in work [9], then the system will be closer to the optimal, the larger the gains of the elements shunted by 
stabilizing devices, 
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ON THE LIMITING SPEED OF ACTION OF SERVO SYSTEMS WITH POWER, MOMENT 
AND RATE LIMITATIONS ON THE EXECUTIVE ELEMENTS* 


E. A. Rozenman 


(Moscow) 


The problem of the shortest transient response is considered for servo systems 
with power, moment and rate limitations on the executive elements. The phase 
plane of the system states is subdivided with respect to the character of the optimal 
transient response for the stated limitations, From the solution of the variation prob- 
lem, the form of the shortest response is determined. 


In works [1-5] the problem of the optimal (shortest) transient response, with limitations on the instanta- 
neous coordinate values, is considered. It is established that an optimal process consists of a definite number of 
intervals, in each of which one of the limited coordinates is maintained at a previously given limiting level. 


The limitation, considered here, on power dissipation (or heating), together with limitations on current 
and rate, leads to a limitation not only on the instantaneous, but also on the root-mean-square value of one of the 
system coordinates. 


1. Posing the Problem 





We shall assume that the energy source does not restrict the possibilities of using executive motors. Then 
the shortest transient response for a given executive element (i.e., for a given inertia of this element) is deter- 
mined by the following limiting characteristics of the motor; 1) maximum permissible motor current ip; 2) 
maximum rate Wm; 3) motor power, i.e., maximum permissible superheating fm. 


We consider a motor, the moment of which is proportional to the current (for example, a dc motor with 
independent excitation), If the induction of the armature circuit is ignored, the state of the executive element 
during the transient response is completely determined by its angular rate w, and angular coordinate g. In 
this case, the problem of the optimal transient response may be formulated as follows: to find a law of current 
variation i(t), compatible with the limitations (i sip, w = Wm, T = Tm), such that the time ty of system 
motion from the initial state (gy = g), w= Wp ) to the given state@ = yp, w = wp) will be a minimum. Here 
7T(0) =T,. 


An analogous problem taking into account only the single limitation on heating (r= T ,,) was solved in 
‘(6]. It was found that the optimal law for current variation was given by the equation 


3 
z= (1—3-)e (1) 
where z, and 6, are parameters of the process which depend on the initial and given (final) states, y is the rel- 


ative time constant of the heating, 8 is relative time,and z {s relative current, If the heating time constant is 
much larger than the electromechanical time constant, the optimal law of current variation is close to a linear 
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one. We shall consider the problem of the optimal transient response when all the limitations enumerated above 
are taken into account, while making the following practically admissible assumptions: 


a) the heat emission of the motor during the transient response is zero; 
b) the static load on the motor is small in comparison with the dynamic load; 
c) the initial state is characterized by the point(g = 0, w = 0). 


For simplicity of computation, it will be assumed that the superheating r starts at its initial value r 9. 
This is equivalent to varying the limitation on the superheating, which now should be taken as the value rg = 
=Tm~To- 


The results obtained in this paper are easily generalized to the case when it is impossible to ignore the 


heat emission of the motor and its load, and also when the initial state of the system differs from an equilibrium 
state, 


With the assumptions given, the motion of the system's executive element is described by the following 
system of equations: 


i os 
——=2_2f, =. @Y, cs = =f. 
a3 a» as (2) 
Here, x, y,Z, # and 6 are, respectively, the relative values of displacement, rate, current, superheating 
and time. 84, is a relative parameter characterizing the heat capacity of the motor. 


It is required to find such a law 2(B8), compatible with the limitations z= zg, y S yg and 9 = 49g, for 
which the time of system motion 68 )y from the initial state (x = 0, y = 0) to the final state (x = xp, y = yp) 
will be a minimum, 


2. The Set of Optimal Transient Responses with Limitations on Power 





We consider the (xp, yp) phase plane of the final (given) system states, The initial state, as was said, co- 
incides with the origin of coordinates. Through each point of this phase plane there passes its optimal trajectory. 
The optimal transient response, with the given assumptions and with account taken of just the one limitation on 
heating, is described by the following system of equations [6]: 


z= 2,(1—7-), y = 2B (1— ge), = ah" ( — 5): (3) 


The parameters of the process, z and Bg, and also the time for the process, 8 yy, are determined from the 
conditions which characterize the final state: 
B 
2B, {1 _ #) = Ye 


+ 5p (1 gee, ) 
2p, (1— M+ =) = Buds. 


The phase trajectory of the optimal process (response) may be determined from (3) and (4) by eliminating 
the variables z and 6 and the parameters 2), 8, and jy. The equation for the phase trajectory of the optimal 
process, 


y (X. y. XR. YF. 8x, Bg) = 0 (5) 
is an algebraic equation of the third degree [6]. 
If the motor’s thermal characteristics (the quantities By; and $,) are given, then the optimal transient 








g=@ 














Fig. 1 


response is unambiguously defined by giving the coordinates (xp, yg) of the final state of the system, i.e., by 
giving a point of the phase plane. Conversely, giving the quantities z, and 8, unambiguously defines the final 
system state (xp, yy) and the time of the transient response 8),. The giving of any one of the quantities which 
define the optimal process divides the phase plane by the corresponding curve. For example, by giving the mag- 
nitude of the relative length (duration) of the transient response, 8 y4 = const, we divide the phase plane by the 
isochronic curves [2] f,(xp, yf. 834) = 0, the geometric locus of points, the most remote from the original state, 
with identical durations of the optimal transient response. However, the form of the transient response for each 
of the isochronic points is different. 


With the condition 84/8 », = const, the phase plane is divided by another curve: 


: 
fe (x. v.55) = 0. (6) 


The transient response has similar forms for all points of this curve. 


In work [6] (Fig. 5) there is presented a family of isochronic regions, constructed on the basis of the equa- 
tion 


- 2 = = 
12xp — 128 yXp¥y + 485 yp — 8), 9,84, = 0 (7) 


(where 6 yy is a parameter), obtained by eliminating z, and 8, from System (4). For different 8 jy, Eq. (7) de- 
scribes a family of ellipses, overlapping and rotated with respect to each other, and with a common center at the 
origin of coordinates, The properties of this family were studied in detail in [6]. 


We now consider Fig. 1,in which are shown the curves, Ba/8 4 = const. They give a family of cubic para- 
bolas 


3 
touching at the origin. The coefficient a depends on the quantity 8 2 yy = q: 


9 (24 — 1) 
t= 7 = ae Ft ? 
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TABLE 1 We point out certain special curves of 


this family. The types of transient response cor- 
responding to these curves are shown in Table 1. 


gz 
es Curve 1 (q = 0) has the equation 


3/9 
4 = de PuPa (10) 


and is characterized by the fact that the initial 
value of the current is z) = 0. The curve divides 
the phase plane into two parts, To the right of 
this curve z, > 0, and to the left z, < 0. Cor- 
respondingly, all phase trajectories of the points 
lying to the right of curve 1 initially set out up- 
wards from the origin, i.e., necessarily pass 
through the first quadrant. Analogously, all the 
phase trajectories of the points lying to the left 
7 A of curve 1 pass through the third quadrant. Curve 
1 is the geometric locus of points of contact of 
the isochrones{cf.[6]), where dy/86yy = 0 and 
8x/88 xy = 0, i.e., points where an infinitesimal 
increase of time does not lead to an extension 
of the isochronic region, 
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Curve 2 (q = 1) has the equation 
/ 9 
p= Ve 8.9, m, (11) 


and is characterized by the fact that during the entire transient response the current retains the same sign and de- 
creases to zero at the end of the response. Correspondingly, at the end of the transient response (3) _ =1,=0, 
“"M 


i.e., the final rate is maximal, Curve 2 is a boundary of the region of monotonic phase trajectories. Its other 
boundary is curve 1. 


Consequently, the region included between curves 1 and 2 is the region of monotonic phase trajectories. 
Curve 2 is the geometric locus of the isochronous points most distant from the yp axis. 


Curve 3 (q = @) has the equation 


3 —— oo 
Y= V 235 %.% (12) 


and is characterized by the fact that the current remains invariant during the entire transient response. The 
phase trajectories of the points lying on curve 3 are second-degree parabolas, Curve 3 is the geometric locus 
of isochronous points most distant from the xp axis. The optimal responses for these points, determined by start- 
ing from the limitations on heating, coincide with the optimal responses which are determined starting from the 


limitations on current. The corresponding isochrones are tangent at these points, Curve 3 lies within the region 
of monotonic phase trajectories. 


Curve 4 (q = $) has the equation yg = 0 (coincides with the axis of abscissas), It is characterized by the 
fact that the transient response is symmetric, and the final state is an equilibrium one, 


If we suppose that xp and yp in Expression (5) may represent the coordinates of any point in the phase plane, 
then this expression will describe the entire set of optimal phase trajectories. To construct them for the problem 
considered here would entail a loss in graphic clarity, since the phase trajectories of different points can intersect, 
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It is possible to visualize indirectly the properties of the families of phase trajectories on the basis of the iso- 
chronous curves 8,/q = const, and of certain other curves, which will be discussed below. 






























3. Subdivision of the Phase Plane According to the Character of the Phase 





Trajectories of the Optimal Processes given Limitations on Rate 





When there are limitations on rate, and for final states sufficiently far from the initial state, there appears 
an interval where y = yg = const. Expressions (3) or (4), defining the form of the optimal transient response when 
there are no rate limitations, are not valid in this case, Therefore, it is of interest to divide the phase plane in- 
to regions where the rate does not attain its limiting value. In order to determine the form of the optimal tran- 
sient response at points in the remaining portions of the phase plane, it is necessary to solve a new variation pro- 
blem wherein the rate limitations are taken into account. 


Thus, to each point of the phase plane individual phase trajectories are drawn from the origin of coordi- 
nates, these trajectories being noncoincident, either in whole or in part. 


We determine that rate, greatest in absolute value, which is attained during the transition from the initial 
to the final state when no limitations on rate exist. The greatest absolute value of the rate Yg may coincide 
either with the absolute value of the maximum rate | ymaxlor that of the final rate | yp|; Yg = | ymax! for 
| ymaxl > | ye] and yg = | yr] for | ymaxl < | yz . The limiting case is the equality 


Yg = | ymaxl =| yel . (13) 


This equality breaks down into two equalities: 


Ymax = YF (14) 
in the case when ymax and yr are of the same sign, and 

Ymax * ~YF (15) 
when ymax and yg are of different signs. 


We now find the geometric locus of points of the phase plane which satisfy Conditions (14) and (15). For 
Y = Ymax, 9y/8 6 = 0, or z = 0. From this, the condition for maximum rate is found: 





B = Ba (16) 
The maximum value of the rate equals 
1 
ymax = > 208 a- (17) 
From (4) the final value of the rate equals 
B 
vp =%By (1 -=M), (18) 
In accordance with (14), we equate Expressions (17) and (18), obtaining the quadratic equation 
B% — 28 By + BR, = 0, (19) 
with the single multiple root 
Ba=By (q=1). (20) 


The optimal transient response for this case is shown in Fig. 2a. The equation of the geometric locus of 
points of the phase plane which satisfy the condition yp = y,nax gives a curve of Family (8) and has already been 
determined [ cf. (11) and curve 2 of Fig. 1). 


Analogously, we find the geometric locus of points of the phase plane which correspond to the condition 
Yr = ~Ymax: From (17) and (18) we find $28, = —z98yy(1—-8)4/28,). Whence 





B% + 2B x48, — BY, = 0. (21) 








Its equation 





AY I MMMM O 

















7 > 


#-£-region 


WY YI JMU as 








The last equation has two roots: Bg 1.3 = 8y4(+ £2 — 1), of which the negative root is not meaningful, 
since for it the rate does not change its sign, which contradicts Condition (15). For this case 


Ba = 0.415 By (q = 0.415). 
The geometric locus of points of the phase plane satisfying Condition (15) is also a curve of Family (8). 


y= — 0.55 VB te, a*p 


is represented by curve 5 of Fig. 1. The transient response for this case is shown on Fig. 2b. 
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Curves 2 and 5 of Fig. 2, where | ymax! = yp, divides the(xp, yg) phase plane into two regions. In the 
region included between these curves, and including the xp axis, the greatest rate is the previously found maxi- 
mum rate, In the remaining portion of the phase plane, the greatest rate is the final rate, characterizing the 
final state. 


Now let the maximum rate be bounded, i.e., 


| yl= ya. (24) 


On the (xf, yp) phase plane we draw the lines y = y, and y = ~yg (Fig. 3), It is obvious that these lines 
isolate the region of states which alone can be attained. We denote this region by the letter H and call it the 
region of admissible states. The portions of the phase plane above the line y = yg and below the line y = ~y, 
are inadmissible, It is also obvious that the phase trajectories of those points of the region for which | yg| > 
>| ymax] , lying below the line y = yg and above the line y = Ya,» are not limited, and that the optimal tran- 
sient responses for these points remain the same as when there were no rate limitations. This follows from the 
fact that for all these points y,= | yr! > ymax. 


At the points where | yg| < | ymax! . on the basis of the position of the points of the final state and the 
straight lines already cited, it is impossible to speak of limitations corresponding to the phase trajectories of 
these points, For this, additional investigation is necessitated. We isolate those points of the (xp, yg) plane the 
phase trajectories of which are not limited by the maximum rate. For this, it is necessary to determine the geo- 
metric locus of the points, the phase trajectories of which have the same maximum rate, equal to the limiting 
value ya. We shall call the corresponding curve iso- maximal to the rate yg. We set 


Ymax = 428 a = ya- (25) 
Solving (25) in conjunction with (4), we find the equation for the iso- maximal curve in the form 


3y2/? — 3y f +1 (26) 
2p 





z= 
where 
1+y4Vi-y 
y (3+ y) 





Y= Ye [Yar = Tp/ Lay Ta = 8ya/ Bua, / = 
(cf. curves AB and CD in Fig. 3). Fig. 3 gives the subdivision of the phase plane into regions in accordance 
with the character of the phase trajectory of the optimum transient response when there are rate limitations. 


In the region ABCD (E-region), the phase trajectories are not rate-limited. The optimal transient responses 
for this region are found by the same rule which is valid for the case when no rate limitations exist. 


In the regions to the right of curve AB and to the left of curve CD the maximum value of rate is the 
greatest, for the limited case. 


In regions DOA and BOC (sub-regions of the E-region), the final rate is the greatest one, Moreover, the 
phase trajectories are monotonic in sub-regions FOA and COE, 


4. Determination of the Form of the Optimal Transient Response when There Are 





Limitations on Heating and Rate 





As was stated, rate limitations alter the character of the optimal phase trajectories only in the regions of 
(Xp. yp) phase space lying to the right of curve AB and to the left of curve CD. For every point of these regions 
the optimal transfer response consists of the following three intervals (Fig. 4): 


1) the interval 8y; = 6, at the end of which the rate will attain its limiting value y,; the path traversed 
during this time equals 


By 
air = \ (1 —B) 2 (B) a8; (27) 





es 





2) the interval Bry = 8, — 8,, during which the rate remains invariant and equal to yg; during this 
time, the path traversed is 
rv = Ya (8,— 8y): (28) 
3) the interval Bf; = 8, — 8, during which the rate attains its given final value yp; during this time, 
the path traversed is 


B 


Zig = Yu Pa— Ba) + | (Po—B)#(P) AB. (29 


6 


The variation problem of the optimal transient response may be formulated in the following manner. To 
find z (8B) such that the time 8yy = 8, during which the system goes from the initial state (0,0) to the final 
state (xp, yp), Will be a minimum with the following limiting conditions: 


B, 8, 
\ dB = ys, \ dB —y —ys, 


8, 8 


\ (Px — B) 2B -+ ya(Bs— Pu) + vs (Be — Bs) +\ (P—B) 2d8 = 2p (20) 
Ba 


0 
8, 


\ ztdB = Bad. 


Let the desired function z(8) during the different intervals of the process be defined in the following way: 


z = 2,(B) for O<CB<(B;, 2 = 2(B) for Pi B< Bs, 
z= 23(B) for Pa<_B< Bs. 


(31) 


Taking into account also that for Zpear = 0 and z,(8) = 0, the limiting conditions (30) may be written in 
the form 


B, Bs 


\ 2,48 = y,, Z,d8 = Yp— Ya 


By 
\ (Bi — B) 2148 + ya + | Prat = Zp» (32) 


8 8 
| #148 + | 2348 = Pao 


The extremal segments 28) and zg8) are found from the condition that the following functional be a 
minimum 


Bs 8, 


vt a ah] zdB + 2s (\ sa + tae) + 


0 (33) 


+ (0 1 —B)sidB + v0 — Pod +f (Bs— 8) 4d), 
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or 
6, 


» = Bs + data (Bn — Ba) + | (azs + Dao + Aa (Pi — B) 24) oB + 
' (34) 


6, 
+ \ Das + Agz2 +, (Bs — B) 25) dB. 


The independent variables on which the variation is carried out are z;, zy, 8,, 8, and By; A,,... Ax 
are the Lagrangian multipliers. Equating the variation to zero, we obtain 


Ay + 2dgz1 +A, (Bi — B) = 0, 
Ae + 2dszs + Ag (Bs — B) = 9, 


Bs 


— AaYa + Aa2s (Bi) + A329 (Bx) + | 2,48 = 0, 


(Az + AsBs — AgBo) Zs (Ba) + As2$ (Ba) = 0, 


Bs 


A + Aaya + date (Bs) + Ase8 (Bs) + Au 2ydB = 0. 
Bs 


(35) 


Equations (35), together with (32), allow one to determine the form of the optimal transient response 
28) and z{8) and all its parameters: A,, Ag, A 3, A4 and 8;, 8, and 8;. Solving the cited system of equations, 
we find 


B 2¥a 
a“ eu(t— 7): “gh Gy (36) 
ta = tq (BF), hae 7 a 





~~ = —_ - 


- ~~ 6a 6 





The durations of the intervals with varying rates equal 


4y;, 
Bit = aaa 33, 9, (1+c*), Bp, = wae a, c(i +c’) (= V 1—%). (37) 


The time of the transient response is defined by the expression 


Ay; 
8, = fas aa9, (1 +o): (38) 
The length of the interval with the constant rate equals 


4y* 


Biv in 98,9, (4 + c*) (? — de — 2). (39) 


In the case considered, the current varies linearly, just as for the absence of rate limitations, in each of 
the individual intervals of the optimal transient response. Its characteristic feature is this, that at the end of 
the first interval and at the beginning of the third, the current is zero, as a result of which the portions of the 
phase trajectories of the optimal response [or the curve y(8)] join smoothly (without breaks), i.e., there are com- 
mon tangents y = yg at these points (Fig. 4). It is easy to show that portions of the extremals 248) and z(8) 
are always parallel for arbitrary values of xp (the larger, certainly, as is made the limitation on rate), 


It also follows from (36) and (37) that the form of the response during the intervals 6 and Bir does not 
depend on the value of xg. For a certain value of xp = Xp,p, portions of the extremals z, and zg join together 
in such a way that the portion z; becomes an extension of the portion z;. In this case, the interval 6,,, disappears, 
and the optimal phase trajectory touches the line y = ya. It is obvious that all points of the (xp, yp) - ae plane 
where such a process occurs are distributed on segments of the curves AB and CD (Fig. 3). 


5. The Greatest Value of Current for an Optimal Transient Response with 





Heating Limitations 





As is obvious from (3), for an optimal transient response, the current varies monotonically, Therefore, its 
greatest value can coincide either with the initial value z,, or with the final value 


(.--e}. 


Zp = Z 


We transform System (4) so that the quantities z and zp will appear as the unknowns. We then obtain 





2y 6x 33,9 
tpt ty = we ty + 2 = 2+ tot +28 = st (40) 
Whence 
1 1 
= (m?n) > + (mn*)?, z ean (m2n)* , (41) 
where 
33.9 2y" 2y? 


631 








It is easily seen from (41) that as the final state approximates to the initial state, z, and zp increase with- 
out bound, 


On the (xp, yg) phase plane we distinguish the set of points, the phase trajectories of which have the great- 
est absolute value of current, equal to zg. This set will consist of points of the four classes for which z, = Za, 


Z, = “Za, Zp = Zag, Zp = “Bae 
On the basis of (41) we obtain 
i 
y= To? [9Buda (242 + Bud.) + 
0 





E V 8162 92 (2yp% + Buda)® — 962842 (2y p20 + Buda)I, (42) 
op= 5 (3 (Aap yy— 388 92) + 
: (43) 





EV 9 (4ae ye, — 36s, 95) 19228 ys @y — Buda]: 


By giving z, and zp the values +z, = const, we distinguish two families of curves on the phase plane (for 
different zg): the family for z, = 4 zg = const and the family for zp = 4 zg = const. To these families should be 
added the isolated curves for z, = 0 and zp = 0, the equations for which were previously obtained. Fig. 5° gives 
the curves forone value of zg. Each of the families cited consists of the two sub-families for which z, = z, or 











Fig. 5 


Sp = Sg, and S, = —Zq Or Zp = —zq, distributed on different sides of the corresponding isolated curves (z, = 0 or 

zp = 0), symmetric with respect to the origin of coordinates, which is the center of symmetry, This follows from 
the fact that Expressions (42) and (43), describing the families considered here, are odd functions of 2, or zp. Each 
of the sub-families consists of overlapping closed curves which, in the limit as z—> contract to a point, namely, 
the origin of coordinates, In addition to this point, which is common to all four sub-families, the following 
characteristic points should be cited : 


*For the sake of graphic clarity, Fig. 5 was drawn up without strict adherence to scale. 








1. The points A and A‘, at which intersect, respectively, curves of the sub-family z, = zg with curves of 
ZF = —Zg, and curves Z») = —Z,g with curves Zp = Za. For Zheat = 0, the points A and A* are distributed on the 
xp axis, symmetrically with respect to the origin, and correspond to the final equilibrium states for which the 
optimal law of current variation is symmetrical. 


2. The points B and B", at which intersect, respectively, curves of the sub-family z, = zg with curves 
Zp = Zg, and curves Z, = —Zg with curves Zp = —zg. The optimal transient responses for these points occur with 
invariant current, i.e., Z) = Zp. Their geometric locus on the phase plane is the well-known curve yp = 


3 
= V 2Bud ax F(9 = 00). 


For a given value | zal = const, there is singled out from the set of curves cited a closed region, 
ABCA'B'C'A (the ®-region), inside which, at any point, the value of current z =z, is exceeded, On the boun- 
dary of this region, at least one value of current (z, or zp), being the greater, equals zg. Inside this region there 
are sub-regions G, G', R and R', each of which is the intersection of two half-regions (shown as cross-hatched 
in Fig. 5). For an optimal transient response, corresponding to any interior point of these sub-regions, there hold 
the inequalities z, > z, and zp > zg. For any point outside these sub-regions but within the @-region, one of 
the values of current (initial or final) exceeds the value zg, and the other (final or initial) is less than zg. Out- 
side of the @-region, the value of current z = zg is not attained for an optimal transient response. 


6. Optimal Transient Response with Limitations on Current and Heating 





We now consider the problem of the optimal transient response with the following limitations; 


|z|< Za, $< 4,. (44) 


By these conditions, the points of the phase plane belonging to the @-region are defined. Depending on 
the coordinates of the final state, the limiting value of current zg may be attained on the initial portion of the 
response, on the final portion, or on both. We consider the case when the limiting value is attained both at the 
beginning and at the end of the response (Fig. 6). To be definite, let z, =z, > 0. During the intervals 6; and 
Bry, the current is defined as; 


z (B) = zg = const for 0= 8 = By, and 2(8) = —zq = const for 6, 3 8 = By. 
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Fig. 6 


Moreover, let z@ ) = 28) for 8,< 8 < 8,. The optimal response will be determined if we can find the 
equation of the extremal section z(8) and the boundaries of the intervals 6,, By and 6. The limiting conditions 
for the case under consideration may be written in the following form: 


Ot hm + | em rp ine 


(more) 











8, (cont'd) 
Za (28,8, + 28,8, — B? — B2 — B2) + \ (B, —B) 2,48 = ate, 


: (45) 
24 (Pi — Ba + Bs) + \ 2*dB ~ Bud. 


The extremal segment 2z,(6) is found from the condition that the following functional be a minimum; 


v = By — AyZa (B: + Bo — Bs) + ha! (2:Bs + 2B2Bs — Pi — Bs — Ps) + (46) 
Gs 
+ Ags (Bi — Ba + Bs) +\ (Ay20 + Xs (Bs — B) Zo + Agee) aB. 
Its first variation is : 
bv = bv (z,) + Sv (B,) + Sv (By) + uv (Bs) = 0. (47) 
Since the corresponding variations are independent, 
bv (zs) = 0, 8v(B,) = 0, 5v(B,) = 0, dv (Bs) = 0. (48) 


Solving (48), together with (45), we find* 


Bs + By ‘ 
#a(B) = % 5a (1 p55) 














(49) 
33y%  3Y 1 
R=> 2 +77 GA 
—_ 9 Pu% 9%p, 5 
8, = 2 22 2 2, +A, (50) 
aw B,%, YF 2 
Bs = 3 2 tats, 
where 
na 
s-y/s 6 + 45°F + 54 ete 4 27 Mt 
“ 
Correspondingly, the durations of the intervals are 
3 B,o 3Yy 
=g-a tzt—e4, 
— —¢ruee A 51 
Br = —6 z 6 E+ A, (51) 
6, 
m= = tae —jA. 





It follows from (49) that, in the intervals Be the current does not attain its limiting value, the current 
varies linearly, the same as in the previously considered cases, 


* The solution is simplified if it is borne in mind that z,(6,) = zg and z,(8,) = —za. 








7. Subdivision of the @-Region According to the Character of the Optimal 





Transient Response with Limitations on the Current 





Suppose that the absolute value of the current may not, or must not, exceed the value zg. The 6-region 
is defined by the set of final states, the phase trajectories of which are current-limited. Depending on the co- 
ordinates of the final state, the character of the transient response differs for different points of the @-region. 
Obviously, the boundary curves subdividing the @-region into shb-regions in Fig. 5 are inoperative in this case, 
since they were constructed without taking limitations on current into account. It follows from (51) that, with 
decreasing Za, the duration of the interval 81 is decreased, but the steepness of the line z,(6) is increased. 
Under certain definite conditions the interval 81; disappears. By substituting 8;; = 0 in (51), we obtain the geo- 


metric locus of points in the phase plane where, with a linearly varying current in an optimal response, the in- 
terval By reduces to zero: 





* 


ype et 2 na Eat g . (52 ) 








Curves (52) for z) = z, and Z) = —zq intersect in points lying on curve 3 of Fig. 1, where the current, dur- 
ing the entire transient response, remains invariant, and equal either to z, or to —zg. They isolate, on the phase 
plane, a closed region (the isochronic region with By = 8y$ af. in the interior of which the optimal transient 
response consists of two intervals during which the current has the values za). Analogously, by substituting 
6; = 0 and By, = 0, we may find the conditions for which the first interval vanishes, 





Zag (53) 


9, 


/ 
ty = oe +) 3 3 Pits 








or the third interval vanishes 








went V SBE ae, . (54) 


It is obvious that the boundaries which subdivide the @-region into sub-regions are just the portions of the 
curves which do not go beyond the limits of the -region, Each of the boundaries defined above divides the @- 
-region into two parts, differentiated by the presence in the optimal transient response of the corresponding in- 
terval, The set of boundary curves (52), (53) and (54) divides the -region into sub-regions with different 
characters of the optimal transient response. It is convenient to use hatching to define which numbers and which 
precise intervals comprise the optimal transient response. 


We shall hatch that side of the curve which faces the region where the interval is present whose condition 
for disappearance is defined by the given curve.. On the basis of the orientation of the hatching it is easy to de- 
termine the character of the optimal response for each sub-region of the @-region. Fig. 7 gives the @-region 
and its subdivisions for the case Zg = 0.5, By; 9g = 1.0. 


Inside the sub-region bounded by the curves 87 = 0, as stated, the transient response consists only of the 
first and third intervals, Moreover, it is characterized by the fact that, for points lying within the sub-region, 
the executive motor cannot be completely utilized for the heating. In the remaining portions of the phase plane, 
the heating utilization is complete, as is attested by the presence of the second interval. In the other sub-regions 
of the @-region, the transient response consists of the corresponding intervals, as designated on Fig. 7. In order 





* The upper sign is valid for z, = 2, and the lower sign for 2) = —za. 
* * The upper sign of + is valid for z) = zg and the lower for z, = —za. 
***The minus sign under the radical corresponds to 2) = Z, and the plus sign to Z) = —Zq. 
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Fig. 8 


to determine the parameters of the optimal transient response for each of its forms, it is necessary to solve a 


variation problem analogous to that studied in Section 6. In all cases, the current varies linearly during the in- 
tervals where |z| < zg. 


Table 2 provides graphs of the optimal transient responses and expressions for determining their parameters 
for the various sub-regions of the @-region (rows 3, 4, 5 and 6 of Table’ 2. 


8. Optimal Transient Responses with Limitations on Heating, Current and Rate 





We now consider the general problem of the optimal transient response in the case when all the limitations 
exist; on heating, current and rate. The solutions of this problem will differ, depending on the mutual orienta- 
tion of the E-region (Section 3) and the @-region (Section 5). It follows from the former that the limiting value 
of the rate is attained for final states sufficiently remote from the initial state. Conversely, the limiting value 
of current is attained for final states sufficiently close to the initial state, 
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The relationships of the quantities By; ,, zg and yg can be such that none of the phase trajectories of 
points in the’@-region will be rate-limited, and none of the phase trajectories of points in the H — E-region 
(Fig. 3) will be current-limited. In this case, the configurations of the E-region and the é-region are mutually 
independent and remain the same as for the individually considered limitations on current and on rate. All pre- 
viously obtained relationships also remain vali¢. ‘ith this, the optimal response cannot, for any point of the 
phase plane, simultaneously contain intervals with current limitations and with rate limitations. 


In order that the E-region and the @-region be mutually independent, in the sense of the previous para- 
graph, it is necessary and sufficient that the @-region be a sub-region of the E-region (cf., for example, Fig. 8). 
A proof of this assertion is given in the Appendix. 


Now, let the line y = y, = const intersect the @-region. In this case, the phase trajectory will be rate- 
limited in some portion of the -region. This portion may be singled out in the same way as this was done in 
Section 3, by constructing the rate iso-maximum ymax = Ya- 


Within the @-region we construct the curve which is the geometric locus of points at which the maximum 
rate equals the final rate, i.e.,points for which ymax = yg. Since the @-region consists of sub-regions with dif- 
ferent characters for the optimal transient response, the curve ymax = yF Will not be a smooth curve (as in Sec- 
tion 3), but will consist of curve segments defined by the relationships which are valid for the corresponding sub- 
regions of the @-region. In the problem considered here, the curve ymax = yf goes through the sub-regions 
I-Iil and I-II (Fig. 7, the dotted line),and consists of three curve segments. In sub-region I-II the condition 


Ymax = yf holds for 6g = 6y. For this sub-region, the equation for the curve segment ymax = yr is written in 
the form 





y 3 B29? 
2+ — dep (55) 
tS ae 2° 


In sub-region I-III] the condition ymax = yg holds for 8, = m, i.e., for z = zg = const. The equation for 
the corresponding curve, 


Ye = 22%, (56) 


is simultaneously the equation for the phase trajectory of point B or point B’, and describes a parabola going 
through these points, the vertices of sub-region I-III. At these points the curve joins with the segments of the 
curve ymax = yF going through sub-region I-II. 
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Let the line y = yg = const, intersecting the @-region, intersect the curve ymax = yr at some point C(C') 
(Fig. 9). This point is simultaneously the initial point of the iso-maximum ymax = Ya (Ymax = ~Ya)- Depend- 


ing on the value of yg, the point C can fall in sub-region I-II or in I-III. Correspondingly, the following cases 
for the path of the iso-maximum ymax = yg may occur; 


1. The iso-maximum goes through the sub-region I-II and reaches the boundary of the @-region at the 
boundary of sub-region I-II (Fig. 9). 


2. The iso- maximum goes through sub-regions I-III, I-II-III and I-11, and reaches the boundary of the - 
-region at the boundary of sub-region I-II 


3. The iso-maximum goes through sub-regions I-III, I-II-III and II-IiI, and reaches the boundary of the 
®-region at the boundary of sub-region II-III. 


To determine the equations of the iso- maximum segments, it is necessary to append to the equations de- 
termining the parameters of the optimal process (response) (cf. Table 2) the equation determining the maximum 
rate. For the sub-region I-II we have 


6 
Ymax = Za (6m en + ys = Ya- (57) 


Solving the system of Table 2, row 3, together with (57), we obtain the expression for defining the iso- 
maximum in parametric form: 











¥s [2a — aut) — 3° 2(u— 1) 
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Here u = 83;/6, is a parameter which must be given within the limits1 < u < 2. 


In sub-region I-III the iso-maximum coincides with the phase trajectory passing through the point of in- 
tersection of the line y = y, with the curve ymax = yr- Its equation is 


Yp = 2y4— 22,Xp . 


To determine the iso- maximum in sub-region I-II-III, to the system of equations on row 5 of Table 2 
should be appended the equation 


8 
¥ max = Za (61 + +) = Ya- (59) 


The corresponding equation for the iso- maximum is obtained in the form 





(60) 


ta (1-98) +209) +3 a 


ss 


where y oh. 3 
a 


Analogously, the equation of the iso- maximum for sub-region II-III is determined by joining to the sys- 
tem of equations of Table 2, row 4, the equation 








Ynax = 7. (Bir — Ba)?- (61) 


The equation for computing the iso- maximum for this sub-region is found in parametric form: 


2, (u—iyP’ (62) 


where u = 8y;/ 6, is a parameter which must be given close to the value u = 2. 
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Fig. 10 


The segment of the iso-maximum from the point where it meets the boundary of the -region to the inter- 
section with the line | y| = y, is determined by the same formulas as in Section 3, The line | y| = yg, together 
with the iso-maxima ymax = Ya 49d Ymax = ~Ya: single out the closed E'-region of the phase plane in which 
the rate is not limited, The optimal transient response for points lying in the E'-region is determined on the 
basis of our previously obtained expressions, In the portion of the é-region which lies within the E*-region, the 
limiting value of the rate is not reached, In the portion of the E’-region not lying in the -region, the limiting 
value of the current is not reached either. 


We now consider the form of the transient response for points of the H — E'-region of admissible states, 
found to the right of the iso- maximum ymax = ya and to the left of the iso-maximum ymax = —yq- It is clear 
that in these regions the limiting value of the rate is reached for any point of the regions. Through the point 
D(D"), the intersection of the iso- maximum with the boundary of the é-region, we draw the line D— D(D'— D") 
parallel to the x axis (Fig. 9). The line divides the region in which the limiting value of the rate is reached 
into two parts. At the points lying below the line D— D, the current does not reach its limiting value. At the 
points lying above line D — D, the current, as well as the rate, is limited. The form of the transient response 
for these points is determined by the fact that through each sub-region of the é-region passes the iso- maximum 
ymax = ya. If from the points of intersection of the iso- maxima with the boundaries of the sub-regions of the 
@-region are drawn lines analogous to the line D — D, then these lines will further subdivide the region where 
the current and rate are limited into sub-regions differing by the character of their transient responses 


For sufficiently small values of y,, the entire E-region may be located inside the é-region. In this case 
the optimal transient response with smoothly varying current cannot occur at any point of the phase plane. 


We now determine the form and the parameters of the optimal transient response for the H — E'-region 
included between the line D — D and the line y = y, (and also for the symmetric region between the line D‘*-D' 
and the line y = ~yg). For this, generally speaking, it is necessary to solve the corresponding variation problem. 
However, by taking into account the results already obtained, the form of the optimal response in these regions 
may be obtained more simply. 


Between the lines D— D and y = yg, we draw some line y = yn intersecting the iso-maximum ymax = Ya 








at the point N. The optimal response for any point lying on this line will differ from the optimal response for 
the point N by the presence of the interval Byy with invariant rate y = y,, where z = 0 (Fig. 10). On both sides 
of this interval, exactly as occurred in Section 3, there are situated 
two intervals, 87; and By, with linearly varying current. The form 
and the parameters of the optimal transient response up to, and 
after, the interval Byy are determined solely by the point of inter- 
section of the line y = yy with the iso-maximum y,yay = Yg- The 
duration of interval 6 yy is determined by the coordinate yr. 








Biv = (xp — x1 — %). (63) 


Ya 
Here x, and x, are, respectively, the variation of coordinate 
x during the intervals of the optimal response before, and after, 
the interval Bry. 








Fig. 11 From what has been said it follows that if the boundaries of 
the E-region and the -region intersect, then each of these regions 
is changed and goes over, respectively, into the E’-region and the @*-region. The E'-region remains closed, 
but its boundaries are different from those of the E-region. The @*-region becomes an open one. To its boun- 
daries, proceeding from the E‘-region, are added the corresponding lines y = 4 yg, D-D and D’-D". The 
possible transient responses with limitations on heating, current and rate are given in Table 2. 


The author wishes to express his gratitude to A. B. Naishulia for his counsel on a number of questions in 
this work. 


APPENDIX 


In order that the E-region and the @-region be mutually independent, it is necessary and sufficient that 
the @-region be contained in the E-region. 











Fig. 12 


The necessity of this condition stems from the following. The boundary of the é-region is defined by the 
absence of limitations on rate, but the boundary of the E-region is defined by the absence of current limitations, 
If the boundaries of the E-region and the @-region intersect, then the boundary of the -region will simultan- 
eously lie in the H— E-region, With this, the portion of the boundary of the @-region leaving the H — E-region 
will also be a portion of the E-region boundary, and will thus change in leaving the @-region, In order to prove 
the sufficiency of the given condition, we first prove the following assertion; if the boundaries of the isochronic 
regions do not intersect, then the phase trajectory of any point inside the isochronic region or on its boundary 
may not go beyond the limits of this isochronic region. 











Indeed, if the phase trajectory of some point M, lying in the isochronic regionI1 (Fig. 11), were to pass 
beyond its boundary, then point N, lying on this phase trajectory but not lying within region I1, could be reached 
in a shorter time than the time of isochronism of I1, which contradicts the definition of isochronism. Consequent- 
ly, if 11 is isochronic with time 8,, and 12 is isochronic with time 6, <8,, then it is always possible to find a 
point on the isochrone I1 whose phase trajectory attains a greatest rate y)4; greater than the greatest rate y)y42 
of the phase trajectory of any point lying on the isochrone 12 (Fig. 12). 


W= now consider the set of phase trajectories of points of the é-region. It follows from what has been said 
that the phase trajectories having the greatest rate may only be attributed to points lying on the boundary of the 
@-region. 


We now explain how the greatest value of rate is reached by the phase trajectories of the points in the @- 
-region. For this, obviously, it suffices to consider only points lying on its boundary, For Zz, = Za = const, the 
greatest value of rate is reached by the phase trajectory for which 8, = 8yy. Analogously, for zp = za = const, 
the greatest value of rate is reached by the phase trajectory for which 8g = 0. From (40) we find the coordinates 
of the corresponding points of the @-region. For point A (Fig. 7) (where 2 = 0), XP, = 385, 89/2), yin = 
= 38 4;9,/2z,; for point B (Fig. 7) (where zp = 0), xs = 384,93 /2z%, Yeas = 38,5, /2z,. 


Substituting the value of yp), in (42) and (43), we convince ourselves that at these same points there is a 
maximum and the boundary itself of the @-region. Thus, the greatest rate, with respect to the phase trajectories, 
of points of the @-region, coincides with the value ypyy of the point of contact of the boundary of the -region 
with the corresponding line ypyy = const. 


Let the limiting value of the rate be y, > yryy- Obviously, in this case the @-region is located inside the 
E-region, since for none of the points of the @-region the limiting value of rate y, is reached. In the limiting 
case, when yg = yp, the boundaries of the E-region and the @-region touch at points A and B (Fig. 8). 


On the basis of what has been said, we conclude that: if the line y = ya = const does not intersect the @- 
~region, then none of the phase trajectories of points lying in the @-region is limited in rate, and none of the 
phase trajectories of points in the H — E-region which reach the state where rate is limited will be limited in cur- 
rent and, consequently, the -region is a sub-region of the E-region. 
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CONVERSION OF OUTPUT INITIAL CONDITIONS IN A LINEAR SYSTEM 
WITH VARIABLE PARAMETERS INTO AN EQUIVALENT INPUT SIGNAL 


A. V. Solodoy 
(Moscow) 


Formulas are deduced to convert the output initial conditions in a linear system 
with variable parameters into an equivalent input signal in the form of a combination 
of various order 6-functions. The obtained formulas are used to determine the output 
transient process, resulting from the initial conditions mentioned above, when the sys~ 
tem's pulse transient function is known. 


1. Introduction 





In the analysis of dynamic systems there often arises the problem of analyzing the excited system for the 
case where the system was in motion prior to the time of application of the input signal. Moreover, at the in- 
stant of application of the input signal, the output function and its n— 1 derivatives have specific values deter- 
mined by the system's motion preceding the input signal's application, These values are the initial conditions 
for determining subsequent motion, and are taken into account in the determination of the system's complete 
output signal. In linear systems with constant coefficients, the initial conditions can be easily taken into ac- 
count if Laplace transforms are used. In this case, the system's transfer function is multiplied by the operator 
equation determined from the initial conditions. This operator equation may be treated as an equivalent input 
signal delivered to a system initially at rest. 


This method is not applicable to linear systems with variable parameters, To take into account the ini- 
tial conditions for such systems it is necessary to determine the arbitrary constants of the general solution of 
the system's equation, which requires knowing the n linearly independent solutions of the homogeneous equation. 
However, this method is not too practical in view of the difficulties associated with finding the fundamental 
system of solutions of the equation with variable coefficients, Moreover, in the analysis of a system it is always 
desirable to have a single approach based on the examination of the action on a system initially at rest, when 
the output process can be determined from the system's basic dynamic characteristic — the pulse transient func- 
tion, With such an approach it is necessary to transform the initial conditions at the output of the system to an 
equivalent input signal, i.e., to transform an excited system to one initially at rest but with a different form of 
input signal, Such a transformation becomes possible if the equivalent input circuit is represented by a linear 
combination of various order delta functions with some coefficients determined by a specific means, 


A linear system with variable parameters is described by the differential equation 
L (P» *) Xout = F (Pp. *) Xin, (1a) 
where xjp is the input signal, Xoyt is the output signal of the system, and 
L(p, th = Ya (t) p’, F (p, t) = 5) 5; (t) P’, p= 4 / at, (1b) 
i—1 j=1 


are linear differential operators. 








Since the right side is a known function of time, Eq. (1a) may be written in the form 
L (P. *) Xout = Yin. (1c) 
where yin = F (p, t) Xjp- 
We use the solution of Eq. (1b) for yin = 5 (t— &) as the system's pulse transient fu:.ction: 
Lip, t)G(t, §)=8(t—&) (2a) 
with initial conditions 
Gi” (t, &) lee. = 0, v=0, 4, 2,..., n—4, (2b) 


where G (t, &) is the pulse transient function and 6 (t— £) is the delta-function effective at the time t = —. 
The delta~function’s fundamental property is determined by the identity [1] 


+c 


#()= \ 1()8e—H at, (3) 


—wo 


while the delta-function itself takes on the values 


nee? (4) 


= * 


In this way, substitution of yin = 5(t— € ), with its meaning taken in the usual sense, will yield a mathe- 
matically unacceptable expression for the right side of Eq. (1c). If, however, this substitution is treated as some 
symbolic expression which, upon integration, yields the transformation defined by expression (3), then the indi- 
cated operation takes on a definite sense. The validity of such a treatment follows, as is well known, from the 
fact that integration and passage to the limit are commutative, which is the basic condition for introducing the 
concept of delta~functions, Let us indicate some of the properties of delta-functions to be used further on. 


On the basis of (3) and (4) we will have 


{()8(¢—§) = 7G) 8 (t 4), (5) 


since, in view of (4), the product f(t)5(t—€ ) is not zero only at the pointt = € , when f(t) = f(&). The ob- 
tained identity is symbolic and its only meaning is that both its sides are equivalent when they are used as the 
right side of Eq. (1c), Thus, multiplying the left and right sides of (2a) by f(E), we obtain L{p, t)Xout = 

= f(E) 6(t— £), where xour = f( E) G(t, &). 


On the other hand, making use of the basic integral relationship for linear systems Xoyt (t) = 
t 


= f G(t, u) yin (u) du and assuming that yj, (t) = A(t) 6 (t— € ),we will obtain 
~ t 


Zour (t) =\GU, u) f (u)8(u—8) du =f (2)G(t, ®). 


0 
In this way, we are assured that the input signal, whether in the form f(€)6(t— €) or in the form 
f(t) 6(t— €), yields the same output signal, which is confirmed by identity (5) in the sense mentioned above. 


Differentiating (5) with respect to‘the variable € , and taking into account the obvious identity 


Se—8) =(— 1 Zae—p, ) 6) 


we will obtain 
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d® a” 
-_——- aw & =i Y mes E — 7 
1(t) Fb e— 8) = (—1) 3 — OL 1 
with identity (7) having the same meaning as identity (5). 


Differentiating both sides of Eq. (2a) with respect to the variable € , and again taking (6) into account, 
we will obtain 


a” 
L , 1G, t, &) = —-4(t — 6), 
aaa (P, )G(t, ) = 3-2) 


G,(t, ) =(—1y E(t, 2 (8) 


denotes the response of the system to an input signal in the form of the v -th derivative, with respect to time, 
of the delta~function. 


Obviously, on the strength of the system's linearity, multiplication of the input signal by any quantity 
which is constant with time t results in the multiplication of the output by the same quantity. 


The rules of transformation for products with delta~functions, defined by identities (5) and (7), must always 
be strictly adhered toto avoid erroneous interpretations during their use. For example, as a result of the incorrect 
interpretation of one of the properties of delta-functions, L. A. Zadeh [2], in the solution of a problem analogous 
to the one being studied in this paper, came to the incorrect conclusion that a varying-parameter system behaves 
itself, in the region of the point t = €, like a system with its parameters “frozen” at the instant € , i.e., like a 
constant-parameter system. As a result, the equations obtained by him are actually only applicable to constant- 
parameter systems, 


2. Conversion of Initial Conditions to an Equivalent Input Signal 





The motion of an initially excited, varying-parameter system for t= € is described by the equation 


L (Pp. t) Xout = Yin 
with initial conditions 


Hut 
dt* 





| , = (5): a=0, 14,....0n—4. (9) 


The input signal yjy acts on the system at the instant t = €,. 


For further investigation, without loss of generality in the results, let us set yin (t) = 0 for all t since, on 
the strength of the superposition principle, the effect of yin(t) on the behavior of the output signal can be ex- 
amined separately, The problem is reduced to obtaining, by means of a specially selected signal,* at the out 
put of a system initially at rest, a process identical to the process at the output of a previously excited system, 
for values of t= —,. The mathematical aspect of the problem consists of finding a right-hand side of the in- 
homogeneous equation 


L (Pp, t) X Sur = Rt, &) (10) 
with initial conditions 
i 
dt* |t—e_ 





=, ami, 2,..., a—14, (11) 


* The point of application of the signal yjn, and not xjn, is considered to be the input of the system. This is 
convenient from methodical considerations, but (if necessary) transition to the point xj, can be achieved quite 
easily, as will be seen from what follows. 
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so that the solution of this equation in x* pyr, being a one-sided function in accordance with the conditions of 
{11), would be identical to the solution of the homogeneous equation 


L (Pp. t) Xout = 0 (12) 
with initial conditions (9) for values of t= £,. 


The one-sided function x*oyt can be easily obtained from the function xoyt by multiplying the latter by 
the unit function 1(t— € ), which will make all values of xoyt , for t< € , equal to zero: 


X* out = Xout 1(t— €). (13) 


Taking derivatives of Eq. (13), up to an including the n-th order, taking into consideration relationship (5), 
and introducing the notation of (9), we will have 


Tout (t) _ Z out () 4 (t git ), 


dz" dz 
—fut = —out 4 (t —§) + x (%)8 (¢ — 8), 
d’z° dz (14) 





sot = ut 4 (t — &) + ay (8) 3(t— 8) + mE) H(t —8), 


d"x° d"z 


Out ae — SUE A (t — &) + tena (8)8 (tf) +> += + (8) 8" (¢— 8). 


Multiplying the first line of (14) by a,(t), the second by a,(t), etc., up to ap(t), adding the obtained ex- 
pressions and grouping terms, we will obtain 
L (Pp, t) X* our — [L (Pp, t) Xoutr] 1(t- E) = RCE, E), (15) 
where L{p, t) is defined by expression (1b), and R* (t, — ) has the following form: 


Re (t, &) = [an (t) ta—a(€) +. . -+ a (t) ay (8) + ay (t) x (E38 (t —8) + 
+ [an (t) tn—a (E) +. « + Gs (t) 4 (8) + ag (t) 2p (E)] 8: (¢ —8) +... (16a) 
vee # [an (t) 2 (€)) Bf"? (¢ — 8). 


For t2=€ the expression for xouyt is given by Eq. (12); hence, L{p, t)x* out = R*(t, €), and the sought~for 
right side of Eq. (10) will be determined by expression (16a), i.e., R(t,— ) = R*(t, €). Transforming this expres- 
sion, with the help of relationship (7), and grouping terms, we will have 


R(t, &) = {{a, (3) — a (&) +... + (—1)” “al (1 2 (&) + 
+ [ag (8) — a5 (8) + ag () —. . «2 (&) +. . + an () tna (8)} 8 (@ — 8) + 1) 
+ {Lag (&) — 2g (&) + Bay (&) —. . .J zy (E) + [as (6) — ag (6) +. . Ja, (§) +... 
«+ et dn (8) Sn—a (&)} 8 (t — 8) +... + (— 1)” “an (&) ty (8) 8)" (¢ — 8). 
Let us introduce notations for the expressions in the square brackets which appear in the multiplier with 
&(t-€) 
Ag (§) = a (&) — a () + a5 (8) —. ..+ (— 1)" 2a!" (6), 
Ay (§) = ay (8) — a5 (&) + 04 (¢) —. . + (— 1)" al (8), (17) 


—— a a - &. Bee & © © a eo eC ew £€ he wo’ se’ O's “o'"> «© 6S 


An—1 (&) = ap (6). 
We write them in operator form: 








Ay (9, &) = a (&) — gas (&) + gas (E) —. . + (— 1)" "9g" tan), 
Ai (q, &) = ag (€) — gas (2) + g*ay (&) —- . + (— 1)" *q” aan), (18) 


SO OS CVE te -- o's "OVO *S G'S" wVet oF Ow ss ae 6¢ ees. = 


An—1 (9, &) = an (8), 
where q =d/dE. 


From a comparison of (16b) and (18) it is easy to find the rule for building up the expression in the square 
brackets, which appear in the multipliers associated with & ()t— ¢),by means of the coefficients in (18). This 
rule consists of differentiating k times, with respect to q, the coefficients of Aj (q, € ) which appear in the multi- 
plier associated with the zero-order delta-function, and then dividing by k! (i.e., the calculation of the expres- 


i at : 
sion = ag A; (q, §)). 
Let us introduce the notation 


Bo (é) =A 0 (9; é) Zo (é) + A, (q, é) ry ” +. ° ‘+ An-1 (q; é) Tn-1 (8), 





Bi (&) = ta (®) + 5 “ie 24) beech = A auaa(t)s 





(19) 
Balt) = 5 (242, 0) + Sn) +...+ ts @] 
sein geveutform OTe ers sree es 
Bx (&) = ‘S =a <* 2, (E)- (20) 
At the same time we will have 
R(t, &) = Bo ()8 (¢— &) — Br (©) & (¢— 8) + Ba BH (¢—8)—... we 


of (— 4)" 82 OR” PO —#. 


Eq. (21) is the sought-after input signal in the system initially at rest, which will bring about in the output 
a process equivalent to the process at the output of an excited system with the initial conditions given in (9). In 
so far as the coefficients 8\(£) depend only on the instant € , they are independent of time t. At the same time 
the output process can be expressed, by means of the system's pulse transient function G(t, t ), on the basis of 
relationship (8), in the following form: 





Zur (t) = Bo (E(t, ) + @) ME! +. + Bana @) EE. (22) 


Eqs. (18), (19), and (21) show that the varying-parameter system does not behave at all like a system with 
“frozen"parameters in the vicinity of the point t = € ,insofaras the coefficients 8, (£) contain the derivatives 
of the coefficients of the equation in aj(t), which would vanish on “freezing.” This is just what happened in 
L. A. Zadeh's work because he assumed / (t) 8{”) (t) = / (0) 8{”’ (t), which is incorrect, 


In the case of a constant-parameter system, the coefficients in (21) will be determined from the equation 


k 
Be =(—1)* S$) ai4s%—« (6), (23) 


i—n—] 











which is easy to see if we note the fact that the derivatives of the coefficients of aj(t) will become zero, The 
Output process will be defined by expression (22), as before. 


Obviously, formulas (20) and (21) can be used to solve the reverse problem, that is, to determine the ini- 
tial conditions which exist at the output of a system initially at rest at the instant of application of an input 
signal of the form given in (21) (more precisely, at the instant t = € ,), which is especially important in modula- 
tion where it is almost impossible to obtain signals in the form of delta-functions. It also follows from the ob- 
tained formulas that for R(t, £) = 6(t—€), i.e., when a signal in the form of a delta-function is delivered to 
the input of a system defined by Eq. (10) and with initial conditions (11), the initial conditions of Eq. (12) take 
on the form 


a(E={ 1 (24) 
a, (é) = n— 1. 


In this way, the pulse transient function G(t, € ) is the solution of the homogeneous Eq. (12) with initial 
conditions (24). 


Example; To illustrate the use of the obtained formulas , let us examine the simple case of the system 
described by the second-order differential equation 


ae () “ou Tay (+) “out 7%)e = (25) 
with initial conditions 


= 2, (E). (26) 


tmE 





dz 
Tout (¢) hog = 2, (8); = 


It is required to find the input signal which is equivalent to the given initial conditions, and to find the 
output process, if the system's pulse transient function G(t, €) is known. 


Making use of formulas (20) and (18), we will have 
Bo (8) = [ai (&) — a, (€)] to (E) + a2 (&) 21 (€), 81 (8) = — a (E) xy (8). (27) 
From this, on the basis of (22) , . 
dG (t, &) 
Zour (ts &) = {la (6) — 2 (8) Bo (5) + aa (E) ws (€)} G (t, €) — ae (E) 2 (€) —2— > ory (28) 


We verify the obtained result by the usual method of finding the arbitrary constants of the general solution. 
It is well known that the general integral of Eq. (25) will have the form 


Xout (¢) = Cypy (t) + Coy(t). (29) 


where gt) and g(t) are the linearly independent solutions of the corresponding homogeneous equation. Differ- 
entiating both sides of (29) with respect to time, assuming t = — , and solving the obtained system for the constants 
C, and Cg, we obtain 


zo(E)%E)— mere) mn )—mE) 9% ©) 
A(é) i= a) 





C= 





From here 
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t — o(t) 0, (E 2 
ee aA 2 2118) 5 + Bm: “©, (30) 





where A(E€) is a Wronskian determinant. 
On the other hand, the pulse transient function can be expressed by means of the formula [3] 
f i (t) a(t) 
as (&) 4 (8) | 1 (8) (6) 1° 


Making use of the well-known relationship from Liouville's formula A‘'/A = —a,/ag, it is easy to obtain 
the expression 


G (t, &) = — 


2 (t) 1 (E) — 1 (t) 2 (E) 
A (&) 


t 4 eer 7 , 
OO - QE 8+ in Be Iee, & (31) 





= a, (6) G(t, §), 








Substituting (31) into (30), we will finally have 


dG (t, & : 


su (t) = {Lar (€) — a (€)} 2 (E) + an (E) 21 (E)} G (t, &) — a5 () m9 6) (32) 


Eqs. (28) and (32) are the same, which is as it should be. 
Received October 27, 1956 
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PLOTTING OF ROOT LOCI OF AUTOMATIC CONTROL SYSTEMS 


N. N. Mikhailoy 


(Moscow) 


The paper deals with analytical, graphical, and semigraphical plotting of root 
loci of automatic control systems. The use of root-finders is suggested to simplify 
the plotting of root loci, The possibility of plotting with the help of automatic de- 
vices is proved and block-diagrams of these devices are given. Other possible appli- 
cations for these automatic plotters are given. 


In the design of automatic control systems today, the so-called root methods of synthesis are often used 
[1-4]. These methods make use of the relationship which exists between the distribution of roots in the character- 
istic equation of a closed system and the exponents of the transient processes, for the synthesis of systems posses- 
sing given control characteristics, In the course of calculations based on such methods, the need arises to re- 
peatedly plot root loci, i.e., the paths of the roots resulting from changes in the amplification factor. Because 
the plotting of root loci is a very time-consuming operation and demands considerable time expenditure on the 
part of the plotter, it is expedient to automate this process, which will lead to wider application of root methods 
of synthesis, 


1. Root Loci of Automatic Control Systems 





As is well known, the characteristic equation of a linearized automatic control system has the form 
KG (p) +1 = 0, (1) 


where the transfer function of the open system KG(p) can be written in two forms: 











Z(p) P™ + by PE + ++ +b prt +++ + hip + by 
(Pp) Q(p) Pp’ +a, Pp” +++ +a, pY+---+ap+a, 
or 
KG(p) = K-22 P=): PT) Pot) og _ 200) @) 


(Pp — 91) (P — 92) - - -(P— 4) ---(P— Gy) Q(p) * 


The real coefficients by and ay have unique values determined from the zeros zi and qi.Generally n = m, 


The roots of the characteristic equation (1) p;,P9..-. » Piss-+» Pn are, in the general case, complex conjugate 
numbers pj = x + jy. 


A change in the amplification factor K of the open system results in a shift of the roots in the complex 
plane, along definite paths, the shapes of which are determined by the function G(p), For example, for K = 0 
the roots of Eq. (1) are situated at the poles of G(p), while for K-> oo they tend to the zeros of this function. 
The geometrical location of all the possible positions of the roots of the characteristic equation, for variations 
in K from 0 to x, is called the root locus of the automatic control system, In so far as K can be greater than zero 
(negative feedback) as well as less than zero (positive feedback), the root locus represents the reflection of the 
real axis of the w-plane on the p-plane by means of the function w = G(p). 
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Fig. 1 


The equation of the locus of a system defined by the function G(p) = Z(p)/ Q(p) has the form 


arg 2) = or, 


where s is any whole number including zero. This equation can also be written in another way 


Z (p) 
1| Q (Pp) | rt 
Here I[Z(p)/ Q(p)] is the imaginary part of G(p) 


(4) 





(5)* 
In the calculation of root loci the most convenient formula is [5] 


n+m Com 
> [sin = >) digua a,)| = 0, 
B=1 Y~1 


(6) 
which is equivalent to Eqs. (4) and (5), Here 6, is the amplitude of vector (p — z), or the amplitude of vector 
Ch 


(p — q) with a negative sign > (IIgtan®,) is the sum of all the combinations of the products of 8 -factors. 
y=1 


The location on the locus of a point corresponding to the given value of K is achieved by means of equation 


Q(p) | _ 1) 
Z(p) | 
or 
Q (Pp) 
e Ese Z (Pp) twl=— = ” 
As an illustration, Fig. 1 shows the root locus of a system described by the function 
Ss p—% aoe p+i1,4 
6)=G=ne—wP—m ~ 





P(p +1)(p+9) * 
The zero z, is denoted by a circle while crosses represent the poles, The arrows indicate the direction of 


motion of the roots as K is varied from —co to +00. The solid line represents the locus of the system with 
* As in original — Publisher, 











negative feedback (K > 0),while for positive feedback (K < 0) the locus is designated by a dotted line. Points 
corresponding to specific values of K are marked on the locus. 
























2. Methods of Plotting Root Loci 





Analytical method. For every concrete system it is possible to write down, with the help of expression (6), 
the equation of the root locus, establishing the interdependence of the coordinates of the points on the locus. 
Choosing the values of one coordinate,the corresponding values of the other are found, a table is drawn up, and 
from it the locus is plotted point by point. The value of the amplification factor K, corresponding to any point 
on the locus, is calculated from Eq. (8). 





For example, let G(p) have one zero z, and two poles q, and qg, i.e., 


io Pp— 2% = p—z 
G(p) (Pp — 91) (P — a2) pP?+ap+a, - 





In accordance with Eq. (6), the suena % x and y of an arbitrary point p on the locus are, in this case, 
related through the expression y[{(x — 2)" + y?-z¥ — ayz, — a,] = 0. 


Hence, the root locus of an automatic control system, for which the function G(p) has one zero and two 
poles, consists of the real axis (y = 0), and the circle with radius r = 2 Z1 + 4%, + a and center at the zero 2, 
(or only the real axis, if z? + a,z, + a) < 0). 





As a second example, let us consider the system with two zeros and three poles 


G (p) and (p — 21) (p — 22) pi p?+bhp+b, 
(Pp — 91) (P — 92) (¢ — 9) p® + asp* + ap + a, © 





Eq. (6), for this case, takes on the form 


y [z* + 26,25 + (3b, + ba, — a,) z* + 2 (bya, — ay) & + bya, — a,b] + 
+ y® (22? +- 2b,2 — b, + ayb, —a,) + y® = 0. 


The root locus, plotted in accordance with this equation for a concrete system (i.e., for given numerical 
values of the poles and zeros), is shown in Fig. 2. 


The building-up of the equation of the root locus for complex systems is already in itself very difficult. 
Ever greater difficulties are encountered in calculating from this equation the tables of coordinates which are 
used to plot the locus. In complex systems (for n + m > 6) it is impossible to express explicitly the dependence 
of y on x, which complicates the tabulation and renders the analytical method of little use for the construction 
of loci. 


Semigraphical method. This method is based on the fact that separate parts of a locus may be determined 
if one makes use of rules which result from well-known [3,4] properties of root loci, The main rules are essen- 
tially the following: 


1, For K = 0 the root locus passes through the poles of G(p) while for K-» op it tends to the zeros of G(p). 





2, Since n— m zeros are removed to infinity, the 2(n — m) branches of the locus also go to infinity, ap- 
proaching the asymptotes forming, with the real axis, the angles 





wl ai om 9 
%2.= (l= 1,2,...,"%—m). (9) 


All the asymptotes intersect at one common point, the center of asymptotes xq, lying on the x axis; 


3 tT 3: , (10) 
im] i=] 
n—m : 





to 
















3. The real axis is always part of the root locus, At the same time, those regions of the real axis to the 
right of which lies an odd number of finite real zeros and roots (including the root or pole situated on the right 
end of the region) correspond to the locus of a system with a negative feedback loop, while corresponding to 
the locus of a system with positive feedback are the regions to the right of which there is an even number of 
zeros and poles, 



















(p + 8) (p + 3.25) 
p(p + 3) (p + 6) 





Fig. 2. Root locus of a system with G(p) = 










4. The points of intersection of the locus’ branches with the real axis (branch points) are determined 
through the real roots of the equation 










(11) 





Z(p)Q*(p) — Z'(p)Q(p) = 0 


where Q'(p) and Z"(p) are the derivatives of Q(p) and Z(p) with respect to p. 






5. The root locus is symmetrical about the real axis. 


That part of the root locus which cannot be found from these rules is plotted point by point. The coordi- 
nates of the points are found graphically by the method of successive approximations, making use of the fact 
that any point on the root locus must satisfy Eq. (4). This means that the sum of the angles of the vectors, drawn 
from the poles and zeros to a proposed point on the locus, must equal sw (the angles of the vectors drawn from 
the poles are taken with the opposite sign). If it turns out that the sum of the angles does not equal sm , then 
another point must be chosen, and so on until the point which satisfies condition (4) is found. The value of the 
amplification factor K, corresponding to any given point, is found as the product of the moduli of the same vec~- 
tors. To speed up the process, the calculation of the sums of the angles and products of the moduli of the vec- 
tors can be carried out on a special rule [3]. 














Let us clarify this method of constructing a locus by taking, as an example, a system with the function 
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In accordance with rule 2, the angles of inclination of the asymptotes, and the coordinate of the center 
of asymptotes, are equal to 
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Fig. 3. Root locus of the system with G = 
: 4 (P)= Sept ip + Sep + 6) 





In accordance with rule 4, the equation for determining the branch points has the form 
(p + 1.1)(4p® + 36p* + 82p + 30) — (p* + 12p* + 41p* + 30p) = 0. 


The branch points are given by the real roots of this equation, pa = —0.722, and pp = —5.52. 


The obtained asymptotes and branch points are entered on the map showing the distribution of the poles 
and zeros of the given function (Fig. 3). Zero z, and poles q,, dg, qs, and q, divide the x-axis up into six regions. 
According to rule 3,the regions ao q4, 43Z,, and q,q,, refer to the locus of a system with a negative feedback 
loop, since to the right of each of these regions is an odd number of finite real zeros and poles (five, three and 
one, respectively). Regions q,q5, Z;4g, and q, « belong to the locus of a system with a positive feedback loop 
(shown with a broken line). 


To find the rest of the locus it is convenient to commence from points lying close to the branch points. 
To verify whether point a belongs to the locus, we connect it with straight lines to all the poles and the zero, 
determine the angles 6 ,,03,0 3,04, and 9, , and calculate their sum (6, +0, +03 +64— 9). In this case the 
sum turns out to be —157,4° # —180°. This means that a is not a point on the locus, and hence we must try 
another point lying further to the left or to the right. When a sufficient number of points on the locus are found, 
they are connected with a smooth line, as shown in Fig. 3. According to rule 5, the part of the locus lying be- 
low the x-axis is plotted like a mirror image of the top part. 


If the value of K, corresponding to some point on the locus, is required, the lengths of all the vectors 
drawn to this point must be measured and multiplied together, with due regard to the scale of the graph. 


ea sins method [6]. The graphical method of finding root loci is based on the construction of so-called 
phase loci G(p). If we designate the argument of G(p) by ¢, and assume it to be equal to ¢;, 
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Fig. 4 





then the geometrical location of all the possible values of p which satisfy the equation 
arg G(p) = #4 (12) 


represents the phase locus of this function for g = ¢. 


Giving ¢ different values, we obtain a family of phase loci for the given function, It is not hard to verify 
that the root locus is one of the curves in this family, corresponding to g =s @ . From this it follows that the 
root locus can be found by plotting a family of phase loci and selecting from them the locus g =s. 


A family of the simplest type of phase loci, satisfying equations 


arg (p — 2;) = $2, arg (=~) = Pq, (13) 


is shown in Fig. 4a and 4b, 











From these two families we can obtain the phase loci of the function G(p) = — corresponding 
1 
to the equation 





i 

arg (p — 2) + arg(———) = 9. (14) 
For this it is sufficient to plot them on the same graph, find the points of intersection for which the sum 

¥z + Yq is equal to g,, and join the points with a smooth curve, as shown in Fig. 4c. By adding on more simple 


types of families we can form a family of phase loci for a system of any degree of complexity, whence it is 
simple to determine the sought-after root locus. 
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Fig. 5 


As an example, Fig. 5 shows the plot of a family of phase loci for a system with the function G(p) = 


= nme . The root locus is designated by the heavy line. 


of root loci by means of root-finders. A root-finder is a device used to determine the roots of 
algebraic equations containing high powers 














Q (p) =anp™ + @n_p"? + +++ ap’ ++. + ap +a, =0. (15) 


The principle of the root-finder's operation is based on the selection of values of p satisfying the initial 
equation (15), For this, the variable p, expressed in polar or rectangular coordinates, is fed to the counting sec- 
tion of the device, which develops a voltage simulating the polynomial Q(p). The zero values of this voltage 
are fixed by the measuring device which indicates the corresponding coordinates of the sought-after roots. 


Depending on the degree of automation in the root-finding process, root-finders are either manual, semi- 
automatic, or automatic [7]. The construction of root loci is most convenient on the latter type of root-finder, 
where the variable p is automatically varied according to a preassigned law chosen in such a way that all the 
points in a particular region of the complex root plane are examined. The distribution of roots for the equation 
being solved is displayed on the screen of a long-persistence cathode-ray tube in the form of points, the co- 
ordinates of which are measured by the operator, Fast-response automatic root-finders permit one to observe 
the motion of the roots as one or another coefficient is variede 


If the function G(p) has no zeros, i.e., 





1 
G = ; 


n—1l 


then the characteristic equation takes on the form 


P” + nap” + «+» +a,p+(a, + K) =0. ay 


Since the root locus represents the geometrical location of the roots of this equation, for variations in K 
between + @ and —o9, then it is clear that by varying the free term (a, + K) of Eq. (17), set up in the root-finder, 
one can obtain the points on the sought-after locus, and even separate parts of it. 


The plotting of root loci by means of a root-finder is carried out considerably more rapidly than by any 
of the other methods examined above. However, it should be noted that existing types of these devices can only 
solve equations of the form of (17), corresponding to the transfer function (16). This makes it difficult to intro- 
duce variations in the location of the poles, and prohibits the determination of loci for systems having transfer 
functions containing zeros. The design of arrangements which do not possess the indicated shortcomings is thus 
extremely desirable. One possible version of such a unit is the one described in paper [7], where the root-finder 
circuit is based on rotary transformers and operational amplifiers, 


3. Automatic Plotting of Root Loci 





Analytical and graphical methods for determining root loci are very cumbersome. But in the design of 
automatic control systems there often arises the need to check, qualitatively, the deformation of the root locus 
with variation of a certain parameter (usually the zero or pole of the open system's transfer function), If the 
methods of plotting loci discussed above were used for this, it would be necessary to choose a series of values of 
the given parameter, to plot a root locus corresponding to each value and, from the obtained set of loci, to de~ 
termine the effect of this parameter. Such an operation would be extremely time-consuming and, moreover, 
it might turn out that the results are unsatisfactory and would, therefore, not be used again. It is apparent that 
the process of constructing loci, which, in itself,is not the object of the design of an automatic control system 
but serves merely as a means of obtaining the system's required characteristics, should take up as little time as 
possible. As a result of this, a device for the automatic plotting of root loci is extremely necessary. 


The possibility of rendering loci plotting automatic is based on the fact that the root locus of an automatic 
control system, having the transfer function KG(p) when in the open state, is a conformal mapping of the real 
axis of the w-plane on the p-plane, through the function w = G(p). If a device was available which permitted 
conformal mapping by means of fraction-rational functions, it could also be used to obtain root loci. 


During the operation of the automatic root-finder, designed at the IAT* AN SSSR [8], it was established 
that, with certain modifications in the circuitry of such a device, it could be used for direct and reverse trans- 
formations from one plane to another, by means of a whole rational function, i.e., for the automatic plotting of 
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root loci of systems characterized by functions G(p) of the type (16). 


It is well known that the formation of the polynomial Q(p) in the counting section of the root-finder is 
carried out in an exponential, trigenometric, or algebraic form. In the first case G(p) is represented by an al- 
ternating voltage of fixed frequency, the amplitude and phase of which correspond to the modulus and argument 
of Q(p), respectively, In the second and third cases the counter produces two voltages, one of which simulates 
the real part, and the other the imaginary part of the polynomial Q(p). From Eq. (5) we determine that for 
the automatic plotting of loci on a root-finder, using a polynomial in exponential form, it is sufficient to fix 
those values of p for which the argument of Q(p) is equal to 0 or . To construct loci on root-finders with poly- 
nomials in trigonometric and algebraic form, it is necessary to fix those values of p for which the imaginary part 
of Q(p) becomes equal to zero. 





[7 } co 




















Fig. 6 


For the last case, the circuit for the automatic plotting of root loci of systems with G(p) functions such as 
in (16) takes on the form shown in Fig. 6. The variable p, which is made to change continuously in accordance 
with the motion of a point along a spiral trajectory in the p-plane, is fed to the counting unit CO, where two 
voltages U and V are produced, representing the real and imaginary part of Q(p), respectively. The voltage V 
is fed to the zero indicator ZI which develops at its output a pulse each time V becomes zero, i.e., when Eq. (5) 
is satisfied. These pulses pass through the noncoincidence gate circuit GC to the modulator of the cathode-ray 
tube CRT. The scan for the CRT is derived from the unit which determines the automatic variation of p, and 
the beam follows the same spiral trajectory along which variations in the p-plane occur. The beam is normally 
blacked out and appears only for as long as there is a pulse on the modulator. As a result, if the spiralling of 
the trajectory of p is close enough, there appears on the face of the CRT a bright line corresponding to that part 
of the root locus which lies in the region of the complex plane being examined. 


The voltage U is fed to the level comparison circuit LCC, which is also fed with a constant voltage pro- 
portional to the amplification factor K. When the real part of the polynomial Q(p), represented by the voltage U, 
becomes equal to K, i.e., Eq.(8) is satisfied, the comparison circuit delivers a pulse to the gate circuit, as a re- 
sult of which a break appears in the locus trace on the CRT, indicating the position on the locus of the point cor- 
responding to the given value of K. By changing the value of the dc voltage fed to the comparison circuit, the 
value of K can be varied between any limits. In this way, besides plotting the root locus itself, the circuit enables 
one to locate points on the locus, corresponding to the values of the amplification factor which are of interest. 
Moreover, by manually varying the coefficients of the polynomial Q(p), formed in the counting unit, one can 
observe on the screen deformations in the locus, and thereby determine the effect of different parameters of the 
system on the root locus. 


Fig. 7 shows the block schematic of a root locus plotter designed on the basis of a root-finder with a poly- 
nomial in exponential form. As in the preceding circuit, the automatically changing variable p is fed to the 
counting circuit and to the scanner for the CRT. All the operations in the building-up of the voltage Up, repre- 
senting the polynomial Q(p), are carried out at the carrier frequency produced by the oscillator G. The measur- 
ing circuit has to fix the times when the argument of Q(p) becomes zero, i.e., when the phase difference between 
the voltages UG and Ug is equal to zero. For this, each of the indicated voltages is fed to its own zero indicator 
which produces a pulse each time the ac signal passes through zero; at the same time the phase shifter PS is in 
the zero state, From the indicators ZI, and ZI, the pulses pass to the coincidence circuit CC which produces a 
pulse if both input pulses arrive at the same time, thereby fixing the zero value of the argument of Q(p). The 
output pulses from the coincidence circuit are fed to the modulator of the CRT and produce on its screen a trace 
of the root locus. Not only root loci can be plotted by means of this circuit. It can also be used to obtain the 
whole family of phase loci of the function in (16). For this purpose the phase shifter PS is used. It delivers to 
the indicator ZI, a voltage from the oscillator G, with its phase shifted by an angle y, as a result of which pulses 
appear at the output of the coincidence circuit only when the following equation is satisfied: 
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arg [(Q(p)] = ¢. (18) 
being the equation of the phase loci, By setting various angles y, by means of the phase shifter, the phase loci 
corresponding to them are displayed on the screen.* The discussed circuits permit one to plot automatically 
the root loci only of systems characterized by transfer functions of the type in (16), i.e., without zeros. To ob- 
tain the loci of systems having zeros in their transfer functions, it is necessary to map the real axis of the w- 
-plane into the p-plane by means of fraction-rational functions of the type in (2) or (3). This problem is con- 
siderably more difficult in view of the necessity of dividing a yoltage representing the numerator of G(p) by 
a voltage representing the denominator, 


These difficulties can be overcome in the case where the polynomials Z(p) and Q(p) are expressed in ex- 
ponential form, if the locus is plotted in accordance with Eq. (5), i.e., fixing those values of p for which the 
difference between the arguments of Z(p) and Q(p) is equal to zero or . The corresponding block diagram is 
shown in Fig. 8; it differs from the circuit in Fig. 7 in that, besides the counting circuit COp, it also contains 
the counting unit COz for forming the voltage Uz. The phase shifter PS permits the plotting of families of 
phase loci, At the same time PS has to change the phase of the transmitted voltage without changing its ampli- 
tude. The zero difference between the voltages Uz and Ug is fixed in the same way as in the preceding cir- 
cuit, at the times when the phase loci equation is satisfied; 


g + arg Z(p) = arg Q(p). (19) 


It should be noted that the phase locus, corresponding to some angle g, is nothing else than the conformal 
map of the straight line passing at an angle g through the origin in the w-plane, into the p-plane, through the 
fraction-rational function G(p). Hence, the described system can carry out the mapping of the indicated straight 
lines, as well as the reverse point transformations. 


The use of a carrier frequency in the root locus plotter simplifies the problem of realizing the locus equa- 
tion but, at the same time, significantly reduces the fast response of the device, It seems, therefore, that quick- 
response locus plotters should be designed on the basis of circuits without a carrier. 


4. Applications of a Root Locus Plotter 


Besides its basic function — the automatic plotting of root loci of automatic control systems — the locus 
plotter can be used to carry out a series of other calculations, 





In the design of automatic control systems, it is often important to know what effect variations in the co- 
efficients ay and by , or the zeros zj and poles qj of the open system's transfer function [9], have on the roots. 
The nature of these effects can be determined from the root locus of the closed system, obtained when the 
parameter in question (ay , by , 24 or qj) is-varied between +c and —oo for a constant amplification factor K. 
The discussed arrangements permit the plotting of such loci. 


* With slight modifications, the circuit in Fig. 6 can also be used to plot phase loci. 








For example, let it be required to find the root loci when the zero z; 1s varied in the system, the transfer 
function of which, for the open conditions, has the form (3). The characteristic equation (1), when solved for 
the incremental value of the varying parameter, has the form 


_ Q(p)+ KZ (p) (20) 
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locus of the system, the traasfer function of which 








is given by the right side of expression (20). To 
ah = {ij (i obtain such a locus it is necessary to introduce a 
few additional components into the circuit of Fig.8. 
The pertinent circuit is shown in Fig. 9. The car- 
Fig. 9 rier signal from oscillator G is fed to the stator of 
4 sine-cosine rotary transformer RT. From one 
rotor winding che voltage goes to COz, from the other toCOo. The tangent of the angle of rotation of the rotor 
of RT represents the amplification factor K. The formation of the voltage KZ(p) +Q (p) is achieved in the sum 
circult £, which is followed by the indicator ZI. Besides KZ(p), the counting circuit COz also produces the 
voltage KZ(p)/(p — 24), i.e., KZ(p) without the factor (p — z4), which is fed to the indicator ZI,. The points 
on the locus are fixed in the same way as in the circuit of Fig. 8. 
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The described device can also be used to plot, automatically, direct and reverse mappings, by means of 
the function w = G(p), of circles centered about the origin, and straight lines passing through the origin. As is 
well known, one such mapping, namely the mapping of the imaginary axis in the p-plane into the w-plane, by 
means of the characteristic polynomial H(p) of the closed system, is the well-known Mikhailov locus used in the de- 
sign of automatic control systems. 


Finally, the root locus plotter may be used for the analysis and synthesis of Fourier series and for the solu- 
tion of transcendental equations, the terms of which can be expanded in power series. 


SUMMARY 


1. The plotting of loci is considerably simplified and speeded up by automatic root~finders, but equip- 
ment of this type which has been available up to now is not equipped to handle the solution of the characteristic 
equations of automatic control systems with transfer functions containing zeros. Since systems having zeros in 
their transfer functions are encountered very frequently, it is necessary to produce root~finders specially for the 
solution of the characteristic equations of automatic control systems in a general form. At the same time, these 
units can be used to automatically plot root loci, 


2. The most effective system of automatic root locus plotting is the one in which the locus (or part of it) 
is plotted as a continuous line, rather than a series of points, on the screen of a cathode-ray tube. This enables 
one to easily examine the effect of various parameters of the system on the shape of the locus. Separate auto- 
matic root-finder units can be utilized in the construction of such root locus plotters. 


3. The proposed block diagrams of root plotters, besides the automatic tracing of root loci, can be used 
to obtain root loci resulting from variations of other (besides the amplification factor) parameters, and to plot 
families of phase loci necessary, for example, for the determination of roots of characteristic equations of auto- 
matic control with lag and distributed parameters [10]. 
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AN ANALOG COMPUTER REALIZING CONFORMAL MAPPING 
FOR AN n-TH ORDER POLYNOMIAL 


V. A. Brik, and S. A. Ginzburg 


An analog computer is described that makes it possible to analyze n-th order poly- 
nomials (n =1,... , 10) with real and complex coefficients. The main purpose of 
the machine is to determine the roots of characteristic equations and plot Mikhailov 
loci. 


The conformal mapping of the complex Z-plane into the W-plane, and vice versa, for the equation 
@p + a,Z + agZ*7 +... 4+ a pZ" = W (1) 


has very important, practical significance in the analysis of automatic control systems, The most useful are the 
mapping of the point W = 0 into the Z-plane (determination of the roots of characteristic equations) and the 
mapping of the imaginary axis of the Z-plane into the W~-plane (plotting of Mikhailov loci). In the course of 
events other mappings prove to be useful too. 


Several electrical analog computers are known which will carry out these transformations. These com- 
puters are classified mainly according tothe method in which complex quantities are expressed in terms of vol- 
tages. At the same time, the following types are distinguished: 1) a complex quantity is represented by a sinu- 
soidal voltage with amplitude and phase corresponding to the modulus and argument of the complex quantity; 
2) a complex quantity is represented by two voltages proportional to the real and imaginary parts of the repre- 
sented quantity. 


On the other hand, similar devices can be separated into those which can be used to transform any curve 
in the Z-plane, and those which can be used to map only certain classes of curves. 


Below appears the description of a machine, built at the TsLEM Mosenergo; which can be used to carry 
out the conformal mapping of polynomials up to, and including, tenth-order ones, This machine can be used 
to map any points or curves from the Z-plane to the W-plane (direct mapping), and certain points (and regions) 
from the W- to the Z-plane (reverse mapping). Complex quantities are represented by fixed-frequency, sinu- 
soidal voltages (50 cycles).Before the polynomial can be fed into the machine ,it has to be subjected to some 
mathematical manipulation. 


Let us express Z in the form 
Z mm pel? = 0g fe? = Pingy® (2) 
where pmax is an arbitrarily chosen positive quantity, 0 <6 <1. Substituting (2) into (1), we obtain (for n=10) 


+ ayp_,, Be? + eee + axe" 


vane eee tH aroh ie. BieetIor = W. (3) 


Dividing the left side of the equation by the largest term (axp Khax)» we obtain the new equation 
Ay + A;Bei* + ApB*ei” + --- + A; Bei” — W, (4) 
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in which the right side differs from the right side of (3) by a constant multiplier only, The coefficients Ay, 
Ay, «+ + » Ayg, Obviously do not exceed unity. 


The polynomial is introduced into the machine 
p in this form. The variables 6 and ¢ are represented 
oe ed by mechanical quantities (angles of rotation of shafts), 


5 
. ¥.) , w : CY) with an unlimited range of variation of the argument 
. 6 - le B er the ra 0=s6=1, 
Sy “e oe gy , while B changes over range 


























4 When 6 varies from 0 to 1 the modulus of the 
rf? variable Z changes from 0 to Pp max, a8 can be seen 
from (2). 


Loew 


Fig. 1. Block diagram of the machine: 1) functional 
converter for the polynomial; 2) phase~sensitive cir- 
cuit; 3) zero element; 4) cathode-ray tube; 5) beam 
deflection control; 6) beam brightness control. The dimensions of the region in the Z-plane 

being examined by the machine can be changed by 
choosing different p »,3x, at the same time recalculating, as shown above, the coefficients A, , Ay, . ~~» Ajg- 
Besides this, the substitution z* = 1/z can be used. After this substitution, an n-th order polynomial can be 
transformed to the form 


The quantity W can be observed on the screen 
of a cathode-ray tube. 


Ag(z*)" + Ay (z* P72 +... + Ay-yZ* + Ane (5) 


The region inside a circle of radius unity and center at the origin in the z*-plane corresponds to the region 
outside a similar circle in the z-plane. Consequently, by analyzing polynomial (5) on the machine (to prepare 
it it is sufficient to change some of the coefficients) all the region in the z-plane previously not embraced is 
analyzed. 


The block diagram of the machine is shown in Fig. 1. Its basic part is the functional converter 1 which 
develops two voltages W and z, The phase~sensitive unit 2 converts the sinusoidal voltage delivered to it, which 
represents a complex number, into two voltages proportional to the real and imaginary components, These vol- 
tages are delivered to the deflection plates of the cathode-ray tube, the screen of which represents a complex 
plane. The electronic zero element 3 operates when its output voltage approaches zero; at the same time the 
beam brightness increases, 


For the upper position of switch S,, the voltage W is fed through unit 2 to the plates of the cathode-ray 
tube, and the locus of point W, in the plane, for changes in the introduced variable z, is displayed on the screen, 
When switch S, is in the lower position the display on the screen is the locus of point z in the z-plane and, at 
the same time, the brightness of the beam increases for W * 0, i.e., the roots of the polynomial are indicated on 
the screen. If an external source (dotted square in Fig. 1) is used to feed into the machine a voltage represent- 
ing the complex quantity reJ9 and, with switch S, in the lower position, this voltage is added to voltage W at the 
input of the zero unit, the bright points on the screen will represent the reverse mapping of the point W = re} 
onto the z~plane, The quantities r and @ are fixed previously. 


Fig. 2 shows the electrical circuit of the machine, For the generation of voltage W in the functional con- 
verter, each term in the polynomial is developed separately, The stators of the selsyns SS are fed by three-phase 
lines, because of which the amplitude of the voltage across the rotor's single-phase winding is constant, while 
the phase is proportional to the angle of rotation of the rotor shaft. By means of a gear train, connecting the 
rotor shaft to the input shaft for the quantity g, the rotor of the selsyn, corresponding to the k-th term of the poly~ 
nomial, turns through an angle ky. The moduli of the coefficients A are set up on the potentiometers A, with 
the switches S taking care of the signs of the coefficients, From the potentiometers A the voltages are fed to 
shaped potentiometers R,, on which the travel of the sliding contacts is proportional to 8, while the shape for the 
k-th member corresponds to the k-th order exponential function (for k = 0 there is no Rs). Consequently, a volt- 
age is developed at the output of Rsk which can be represented by the complex quantity A,B*e/(%r+**), The angle 
@o is determined from the mutual positions of the selsyn's rotor and stator when y, = 0. Usually g, = 0, but if a 
coefficient Aj, turns out to be not real but complex then, by turning the selsyn stator (i.e., by changing ¢,), the 
argument of the complex coefficient can be taken into account. We shall henceforth assume that all the coef- 
ficients are real, 
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Fig. 3. General view of the machine: 1) selsyn unit; 2) rheo- 
stat unit; 3) electronic unit. 


The sum of all the voltages, corresponding to the terms of the polynomial, forms the resultant voltage W. 


The functional converter, together with the first-term (A,z) selsyn SS,, and the linear rheostat Rp, are also 
used to generate the voltage +z = +Be/% (the sign depends on the sign of coefficient A,). The introduction of 
the quantities 8 and ¢ into the functional converter is achieved either manually, by means of the handles HA, 
or automatically, by means of motors MM and MA. Contacts Ks and Kg are used to switch in the motors, and 
the direction of rotation is altered by means of switches S, and S,, while the speed is controlled by autotrans- 
formers. 


The shaped potentiometers R, (as well as potentiometers A) are the ring type. 


Exact readings of the variables 8 and g are displayed on the modulus scale and argument scale which 
are connected to the 8 and g input shafts. In this wav, the variable z can be observed either on the screen ofthe 
cathode -ray tube (coarsely) or on the scales (accurately), 


The functional converter (see general view of machine in Fig. 3) is built in the form of two units; the 
selsyn unit which contains the selsyns together with the functional transmissions, motor MA, and potentiometers 
Ag, Ay... ~ » Ayo, and the potentiometer unit containing the shaped potentiometers Rpy, Rng, « - - » Rnio» Rp» 
together with the functional transmissions and motor MM. 


The shaped potentiometers consist of PEK-0.2 type wire wound on notched, ring-shaped forms; they thus 
yield piecewise linear approximations of an exponential function, Each segment of the broken line is adjusted 
by means of a shunt connected in parallel with the corresponding step, so that the total error of the approxima- 
tion does not exceed 4 0.25%. 


For those gear trains in the functional converter which influence the accuracy of the whole unit, back- 
lash is eliminated by utilizing sectional backlash-free gears. 


The electronic unit (Figs, 2 and 3) contains the phase-sensitive circuit, the zero element, and the cathode~ 
ray tube. 


The phase-sensitive circuit utilizes 6N8 tubes — tubes V, and V, in Fig. 2, After a preliminary amplifica- 
tion in the stage consisting of tube V,, the voltage from the output transformer Tr, of this stage is fed to two 
push-pull stages (tubes V,, Vs, and V¢, Vz) controlled by two “electronic keys," made up of tubes V, and V,. 
During the positive half-cycle of the voltage, taken off Rgs, tube V, conducts and the current through it produces 
a drop in the voltage across Rgg, Which cuts off V, and Vs. During the negative half-cycle V, is cut off and V, 
and Vs conduct. V4, Vg, and V, operate in an analogous manner, The controlling voltages across Ca, and Reg 
are 90° out of phase. Transformer Tr, is fed from the ac mains. As a result, the constant component of the 
voltage difference between the plates of V, and V, and of V, and V, is proportional, respectively, to €¢¢ cos a 
and e, sin a, where é¢ is the amplitude of the voltage across the input of the phase~-sensitive circuit, a is the 
phase difference between this voltage and the reference voltage across Rgg. The reference voltage is in phase 
with the voltage on the winding of the stationary rotor of selsyn SS, which is used to develop the constant term 
A, of the polynomial. In other words, the constant components of these two output voltages represent the real 
and imaginary parts of a complex quantity. The indicated constant components are shaped in two RC~filters 
and fed to the plates of the cathode-ray tube. 


The zero element is built around tube V,(6N8). The output of this zero element is connected to a relay 
R, in which contact K changes the brightness of the beam, 





The beam Jeflection (zero setting) is achieved by means of potentiometers Rg and Ry, Ry and Rg form the 
focusing contrci, and the brightness is varied by means of potentiometer Ryy. Ss and Rg are used to change the 
sensitivity of the zero element. The signal fed to the phase-sensitive circuit is regulated by means of S, and 
Rg Or Ry. Also, the separate control of amplification of real and imaginary components is achieved by means 
of coupled potentiometers connected to the circuits of tubes V_, Vs and V¢, Vz. 


The upper position of switch S, in Fig. 1 corresponds to the upper position of S$; in Fig. 2, while the lower 
position of S, in Fig. 1 corresponds to the second and third positions of S; in Fig. 2. The second position of S, 
in Fig. 2 is used during observation of the z-plane for A, > 0, while the third position is for Ay < 0. The cathode- 
ray tube (13L036), with a 130-mm screen diameter, is of image-persistence type (for a properly darkened screen 
the trace is visible for 3 to 4 minutes). This permits one to observe simultaneously all the loci of the variables 
W orz. The tube’s electrodes are fed from the power supply comprising tube V4, (2Ts2S) and transformer Trg. 


The functional converter and electronic unit possess the following basic parameters; 10 v maximum out- 
put from each formed potentiometer, 110 vy maximum summed output voltage, 0.7 v minimum input voltage 
to the phase~sensitive circuit, to produce full deflection of the beam. 


If switch S, (Fig. 1) is in the upper position, the machine displays on the screen the direct mapping of that 
curve along which the operator changes the independent variable z. If S, is in the second position, then the re- 
verse transformation of the point W = 0 onto the z-plane appears on the screen as bright flashes. If the zero 
element is spoiled, then, instead of the reverse mapping of the point W = 0, the region in the W-plane, lying in- 
side a circle of some radius and center at the origin, will be mapped. 
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Fig. 4. Evolute of the variable z and the determination of roots by means 
of the zero element. 








In this way, the plotting of direct transformations is carried out by the machine directly from the intro- 
duced values of z, the transformations of which are required. However, a reverse transformation is plotted by 
the method of plane scanning, i.e., the variable z is changed by some means until the required magnitude of 
W, the reverse transformation of which is required, appears at the output. 


Fundamental practical problems which are solved on this machine are the determination of the roots of 
a polynomial and the plotting of Mikhailov loci. For the determination of a Mikhailov locus a direct mapping 
is carried out; the condition y = 90° is set up, 8 is varied between 0 and 1, and the required W-curve is displayed 
on the screen, 








The problem of determining roots, when z is introduced automatically, is solved in the following way. 
By means of motors MA and MM the variable z changes in such a way that the point z explores all the region 
being examined. For this 8 changes repeatedly and quickly from 0 to 1, while g changes slowly from 0 to 360°. 
The result is a radial evolute embracing region z (Fig. 4). As z approaches one of the roots the zero element 
operates (the operating voltage is equal to oa= | W| oper) and a point is fixed on the screen, During a prelimi- 
nary investigation the zero element can be rendered less sensitive and a rarer evolute can be used to find the ap- 


proximate location of the roots. After this, increasing the sensitivity and utilizing a denser evolute, the values 
of the roots can be determined more accurately. 


A convenient method of finding roots is the one in which the variables 8 and gy are changed manually 
by means of the handles HA (Fig. 2). In this way the operator “feels” for every root. The method of root find- 
ing is, in this case, based on the argument, which arises as a consequence of Cauchy's theorem, that when z 


moves along some closed contour, W rotates about the origin as many times as there are roots within the closed 
contour in the z-plane. 


The manual location of roots is carried out in the following way. A count is made of the number of times 

W rotates around the origin when z moves along the circle 8 = 1, This number of revolutions will indicate the 
number of roots lying in the section of the Z-plane 
being investigated. Next, by making the point z 
travel along slowly widening concentric circles, 
starting at small values of 8 and increasing up to 

@® B = 1, each of the roots is localized. Every time a 
root is traversed the beam (point W) passes through 
the origin. 








¥ Multiple roots are found in the same way, i.e., 
by counting the number of times W rotates about the 

origin when z travels around the root, the multipli- 

city of which must be verified. Upon locating each 

root inside the circle 6 = 1 the operator changes 

Jf P max or introduces the substitution 1/z and locates 

the remaining roots, if they exist. 








Fig. 5. Mapping of the outline of the left z-half-plane If only the system's stability has to be verified 
onto the W-plane. then, instead of finding each root, it is sufficient to 

plot a Mikhailov locus or to count the number of 
revolutions of W corresponding to z travelling around the left half-plane. Fig. 5 shows, for example, the plot, 
obtained on this machine, of the mapping onto the W-plane of the outline of the left z-half-plane for a fourth- 
degree equation. Four rotations about the origin indicate that all four roots lie in the left half-plane and hence 
the system is stable. 


An experienced operator can find all the roots, lying within a given pray in 1-5 minutes. 


The time required to plot a Mikhailov locus is equal to the time necessary to change 8 from 0 to 1, If 6 is 
introduced automatically this time is 1 to 5 seconds, if it is introduced manually the time is 7 seconds or more, 


The accuracy with which problems are solved on the machine depends, firstly, on the accuracy with which 
W is developed and, secondly, on the equation itself. In particular, the accuracy with which the roots of an equa- 
tion are determined depends to a great extent on the location of the roots in the z-plane; the closer the roots 
lie to each other, the less accurately can they be determined. 


From the circuit in Fig. 2 it can be seen that voltage W is the algebraic sum of eleven output voltages from 
identical circuits, representing the members of the polynomial. It can be assumed that each of the summands is 
developed with an error not greater than 2% (the error results from errors in selsyns, potentiometers, and other 
components). This permits us to guarantee the following accuracy of W: the length of the error vector 6 W (i.e., 
the difference between the W built up in the machine and that which actually corresponds to a given z) does not 
exceed 2% of the sum of the moduli of the polynomial's terms for a given value of z. This sum is equal to 

10 


W’ = > | A; | B‘. Therefore, to determine in practice the accuracy with which a particular value of W is developed 
i=0 
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in the machine for a given polynomial and a given value of z, it is necessary, without changing 8, to establish 
the equality gy = 0, to give all the coefficients the same sign, to measure directly on the screen the magnitude 
W* and to calculate 2% of it. This will be the maximum possible absolute error in setting up the polynomial 
&Wmax. This assumption was justified by solving many standard problems. From this follows a method for de- 
termining the accuracy of a calculated root, i.e., the region around the root zy, found on the machine, within 
which the true value of the root is definitely confined. 











-Q 


Fig. 6. Mikhailov locus. Solid line — calculated; dashed — experimental. 


It is necessary to measure W" for the point zk, calculate 2% of it and then, manually varying 6 and ¢, 
locate the contour of that region around the point z, for which the magnitude of the modulus of W does not exceed 
the determined value of Wmax. This zone can be located by means of the zero unit if its sensitivity is adjusted 
so that it operates for a voltage | W| = 0a = 5Wmax (Fig. 4). In other words, the true root lies within a region 
which is the reverse mapping of the interior of a circle with radius 6 Wax and center at the origin in the W-plane. 
It is obvious that while the magnitude of the radius 6 Wmax depends mainly on the properties of the machine, the 
nature of the reverse transformation of the indicated circle from the W-plane is determined by the type of poly- 
nomial. The smaller dW dz| 2=Zk is, the larger is the region in the z-plane which is the reverse mapping of the 
indicated circumference and, hence, the smaller is the accuracy of the determined root. An especially large error 
arises in the location of a multiple root, since at that point dW/dz = 0. Moreover, at points where dW/dz = 0 
the mapping is not conformal. 


If the multiplicity of the root is k, then all the derivatives from the first to the (k — 1)-th become zero. 


The higher the multiplicity of the roots, the greater is the error in determining them. The determination of 
multiple roots can be made considerably more accurate by transferring the origin into the vicinity of the multiple 
roots, by replacing the variable z, in the polynomial, by z' = z— u where, in the general case, u is a complex 
quantity lying near the multiple roots. As a result of these transformations we obtain a polynomial in z', with com- 
plex coefficients, which can be analyzed on the machine. By choosing the corresponding p max it is thus possible 
to separate these roots far from each other and, therefore, to determine them more accurately, In certain particular 
cases u can be made a real number, which leads to a polynomial in z* with real coefficients, which can be ana- 
lyzed 6 on the machine in the usual manner. This is convenient when the polynomial contains conjugate roots, the 
real part of which is much greater than the imaginary, or near-real roots. 
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The examples presented below give a good indication of the accuracy with which problems are solved on 


1. Determination of roots of equation 








Value of root 
N Exact value of root chanbuidl with sendin 
{ 0.200 0.188 
2 0.424  ¢4135.0° 0.425  ¢/135.8° 
3 0.424 ¢525.0 0.425  ¢22A.0 
4 0.671 ¢/88-5° 0.675  ¢63-4° 
5 0.671  ¢/2%6.5° 0.676 1296.5" 
6 0.825 #43 0.835  ¢/13-5° 
7 0.825 2945.7" 0.827 1346.0 
. 0.900  ¢/180.0° 0.902  ¢f!79.8° 
9 0.971 ¢f0t-8° 0.974  ¢fi01.5° 
10 0.971  ¢/%8.2° 0.969  ¢258.0° 








2. Plotting of Mikhailov locus for the polynomial 


0.15 + 0.52427 + 0.85 2° +z*=w. 


Fig. 6 shows the calculated locus and the one obtained on the machine. The differences between the 
calculated and experimental curves are due mainly to errors in the electronic unit, 


The described machine was built at TsLEM Mosenergo*in duplicate and is being used in the dynamic 
model laboratories of MEI** and VNIIE MES*** for the solution of problems connected with the stability 
of power systems. 


Received July 4, 1957 
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AUTOMATIC ZERO DRIFT COMPENSATION IN ELECTROMETRIC AMPLIFIERS 


D. E. Polonnikoyv 


(Moscow) 


The results of investigation into the automatic compensation of zero drift 
in electrometric amplifiers are given. The sources of error with such compensation 
are shown, and methods of eliminating them are cited. Circuits for drift compen- 
sation are presented—circuits which provide minimum error. A description is given 
of an amplifier, developed at the IAT “AN SSSR, with dynamic condensers and with 
automatic zero drift compensation. 


INTRODUCTION 


The problem of measuring small, slowly varying currents (10~"* to 10°** amp ) and emf sources with large 
internal impedances is particularly important in the measurement of large resistances (greater than 10"° ohms), 
ion concentrations, semiconductors and insulators, in mass spectroscopy, ionization and bioelectrical measure- 
ments. 


The measurement of small currents turns out to be necessary, not only in the execution of scientific in- 
vestigations, but also in industry, particularly in relation to automation and to the use of transducers employing 
radioactive isotopes. 


The very best instruments, allowing measurements to be made of small currents and emf sources with high 
internal impedances, requiring very careful technological preparation, are far from satisfying contemporary re- 
quirements. For example, the “Cary model 31" [1] electrometer, developed in the USA in 1956, although pos- 
sessing a sensitivity close to the theoretical limit, still possesses a zero drift of the order of 200 microvolts a day, 
which, with an input impedance of 10 ohms, corresponds to a current of 2° 107"* amp.** Therefore, before 
starting to make measurements, one must establish zero and, consequently, the instrument is not useful for con- 
tinuous control over a long interval of time. 


The reduction of zero drift is a fundamental problem in the measurement of small currents. 


At the present time, the least value of drift, and the greatest current sensitivity, are obtained in electro- 
metric amplifiers with dynamic condensers [2,3], The principal cause of zero drift in dynamic condensers is the 
lack of a constant contact potential between the active electrodes. Even gold-plating the dynamic condensers, 
or immersing them in inert gasses, does not allow the zero drift to be reduced any lower than 4100 microvolts 
per day. However, there is another way to solve the problem, which consists in using automatic zero drift com- 
pensation. The idea of compensating zero drift was presented back in 1947 by Prinz [4], who succeeded in re- 
ducing the drift in a dc amplifier with ordinary vacuum tubes to 1 millivolt, Today, zero drift compensation is 
widely used [5-7] in dc amplifiers with low and medium levels of input impedance (up to several megohms). 


We shall analyze designs for the automatic compensation of zero drift by the use of contact switches, ex- 
plaining what errors are introduced by the switches, and giving methods for decreasing these errors, Based on the 


* Institute of Automation and Telernechanics. 
“ Increasing the input impedance theoretically permits the zero drift to be reduced, but the use of impedances 
greater than 10" ohms is unsatisfactory, since these possess lower stability and introduce additional errors. 
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investigation described here, two types of electrometric amplifiers were developed, the drifts of which approxi- 
mated to the noise level, and the sensitivity of which approached the theoretically attainable limits, which are 
governed by the thermal noise in the input circuit. 


1. The Errors Arising with Automatic Compensation of Zero Drift, 





and Methods of Reducing Them 





One of the simplest schemes for automatic compensation of zero drift is shown in Fig. la. Under the 
stimulus of pulse generator PG, relay P periodically operates, and its contacts Pa and Ph are switched from 
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Fig. 1. Circuits for zero drift compensation: a) with partial compensation; 
b) with complete compensation. 


state 2 to state 1. In state 1 the dc amplifier with gain K is shunted by a 100 % negative feedback loop, and 
the condenser C, is connected in parallel with the amplifier output. With a zero shift Uzs present at the input, 
condenser C, is charged to the voltage Uz, [1— 1/0 + K)]. If we assume that condenser C, is self-discharging, 
that the variation in zero shift is negligibly small during the time interval between generator pulses, and that 
the contacts themselves do not introduce parasitics, then the output voltage in contact state 2 is defined in the 
following way: 





” Uzs ) 
Uout (vin + i+K K. (1) 


It follows from this that the zero shift introduced at the input is unboundedly decreased with increasing K. 
The zero shift may theoretically be reduced to zero even with K not very large if the circuit is modified in the 
K — 





manner shown in Fig. 1b. It is easy to show that for B = the “memory” condenser C, will be charged 


exactly to the value Uzgs, i.e., will implement complete compensation. 


In fact, the zero drift cannot be made arbitrarily small, due to the parasitic signal U, generated by the 
contacts in switching. A special investigation [8] enables one to explain that Uc consists of three components. 
The first of these is engendered by the variation in mechanical tension on the insulators supporting the contacts 
of Ph, due to which leakage occurs on their surfaces. The size of this leakage can be as much as 10°* coulomb 
and is determined by the insulator material and the pressure variation during switching. Obviously, to decrease 
this component, it is necessary to decrease by all means the difference in pressure on the insulators in contact 
states 1 and 2, and to choose high-quality insulator material. 


The second component is related to the transfer of the contacts in an electrostatic field which, due to the 
presence of random charges on the insulator and contact surfaces, generates a potential difference between the 
inhomogeneous surfaces. To decrease this latter component of error, it is necessary to minimize the distance 
of travel of the contacts, and to decrease the electrostatic field strength in the region of contact motion. This 
latter is attained by careful insulator shielding using metals with minimum contact emf with respect to one 
another, and by increasing the distance between the contacts and other parts. 


The third component stems from the presence of contact emf differences between the active surfaces of 
the contacts, leading to the appearance of a charge on one contact with respect to the other after contact break- 
ing. We have established that the magnitude of the third component depends only slightly on the type of contact 








TABLE 1 








Alloy Ag 99:9 Pd 60, Ag 40 Au 98, Cr 2 Au 04, NI5, Mn 1 
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material (cf, Table 1), but quite strongly on the purity and homogeneity of their surfaces. The results given in 
Table 1 were obtained with contact point specimens prepared from 1 mm-diameter wires that had been washed 
with alcohol. All the specimens were tested under strictly uniform conditions by means of a high-sensitivity 
electrometer. The very best results were obtained with contacts of a gold, nickel and manganese alloy which, 
with careful cleaning of the surface, allowed the generated charge to be reduced to a magnitude of 3 « 107** 
coulomb i.e., to 10 microvolts with a switch input of 30 micromicrofarads capacity, With this, the total para~ 
sitic signal was of the order of 15 microvolts. The variation in magnitude of this signal observed during a 
period of several days did not exceed 5 microvolts (1.5 * 107'* coulomb ), or 2.5 * 10™ amp for a commuta- 
tion frequency of one per minute, 


Thus, with a proper choice of material and design,the errors introduced by contact switching can be re- 
duced to a very small quantity, allowing their use for zero drift compensation even in the most highly sensitive 
amplifiers with dynamic condensers. 


We remark that, at present, there is no other means of compensation known which possesses sufficient 
stability to allow zero drift compensation in electrometric amplifiers to be carried out. 


The circuit of Fig. 1, due to Prinz, possesses essential drawbacks which lead to the appearance of addition- 
al errors, The principal drawback is that, with switching of the contacts, the voltage at the amplifier input 
periodically varies from Ujp to Uggs, and overcharging of the input capacitor Cjp occurs. Such overcharging 
leads, on the one hand, to a decrease in sensitivity, since additional energy is required from the signal source; 
on the other hand, the duration of the transient response increases, which requires a significant increase in the 
time Tc for drift compensation, i.e., a lengthier switching off of the input signal. This entails an increase 
in error, engendered by the periodicity of the measurements. In order to decrease these errors, it is necessary 
to minimize the compensation time, 


Moreover, switching off of the input signal for compensation leads, in the circuit of Fig. 1, to undesirable 
transient responses in the measuring or registering instrument connected to the amplifier output, 


Errors also arise due to the possible variations in the zero shift Uzs, in the periods between compensating 
action. However, if the insignificant speed of variation of Uzs is allowed for, and if the compensation period 
Tc is chosen to be not too long, then this error can be made less than the noise level, 


We note that another reason for choosing Tc not too large is the attempt to decrease the error generated 
by the self-discharging of condenser C,. The zero shift can attain several millivolts, but the admissible error 
desired is of the order of a microvolt; therefore, the time constant for the self-discharge of the condenser must 
be approximately a thousand times larger than Tc. To compensate for drift once a minute (Tg = 60 sec ), the 
condenser must have an isolating resistance of not less than 6 * 10" ohms/ microfarad, which is attained only 
with polystyrene condensers. 


We have developed several variations of compensation circuits in which are attained significant decreases 
in the errors connected with the switching of the amplifier to the drift compensator. Fig. 2 shows the block 
schematic of an amplifier which provides complete, for practical purposes, continuity of measurement, The 
device consists of two identical electrometric amplifiers I and Il with gain of K, (each), an output stage III 
with gain K,, three memory condensers C,, C, and Cs, a carrier frequency generator CG for supplying the 
dynamic condensers and demodulators of amplifiers I and I, pulse generator PG, periodically exciting relay P. 
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Fig. 2. Circuit for zero drift compensation which provides continuous 
measurement of the input signal. 


The circuit works in the following manner. In the half-period when contacts are in the upper position amplifier 
I amplifies the input signal and operates on output stage III. At the same time amplifier II is shunted by a 100% 
feedback loop by which condenser C, is charged up to a voltage equal to Uzs by amplifier II with an accuracy 
Uzs 
fitK, 
and Pf are raised), the role of the amplifiers are exchanged: II amplifies the input signal and I is switched to 
drift compensation. With this, Cis charged to Uzs by amplifier I,while C, is connected in series to the input of 
amplifier If and compensates its zero shift. C, is used to retain the signal during the transfer of the contacts. 
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. In the succeeding half-period, after switching of the contacts (Pg, Pp and P, are lowered, Pq, Pe 
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Fig. 3. Circuit for zero drift compensation with selected moments 
of compensation, guaranteeing minimum error. 


Obviously, such a design for an amplifier practically completely eliminates the errors of the circuit of Fig. 1 
which were generated by the periodic switching off < the input signal. In the circuit of Fig. 2 the potential 


iD » where K,K, on >1, is negligibly 
2 


+ 





differences between points a andb, equal to Rs 
fos 1+ KK Tet R 
small, Therefore, when contact switching occurs, practically no overcharging of Cin occurs and, consequently, 
there is no lengthening of the transient response and no decrease in sensitivity. Among the disadvantages of this 
circuit are the large number of contacts and the necessity of having two amplifiers, so that its use is expeditious 


only in those cases when the input signal changes quickly and must therefore be continuously observed. 


Significantly simpler is the circuit of Fig. 3, in which are employed only one electrometric amplifier and 
half the number of contacts, The circuit is analogous to one channel of the device considered above. The 
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distinguishing feature of the circuit is the control of the length Tc of the generator PG cycle by means of 
the introduction of a signal proportional to | dUjp/dt| . This control is so effected that moments of time for 
switching the amplifier to drift compensation are chosen when the speed of variation of the input signal 
dui, | ~ 1 4 ow 

dt Ki Kg dt 
long time, then compensation will also occur, this being necessary because long functioning without zero drift 
compensation would lead to larger errors from this cause than would a short-term switching off of the input sig- 
nal. The circuit for controlling the generator could be implemented as shown in Fig. 4. It contains a two half- 
period rectifier, the output of which is the modulus of K,K,Uin, 4 differentiating network R,C, an amplifying 
stage using tube L, and a generator using tube L, with relay P. 











is a minimum. If the input signal varies continuously with significant speed for a 
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Fig. 4. Circuit for controlling the time interval between moments of 
drift compensation as a function of the magnitude of | dUjn/dt| . 


dUin 

dt 
(with time constant C,R,), the grid potential is increased until relay P operates. With this, the current through 
the tube (L,) rapidly increases, since resistor Rg is shorted by contact Pg. Then the grid potential is reduced due 
to the overcharging of C, via Rs, with time constant CRs, until the relay releases, etc, Contacts Pg, Ph and Pc 
(Fig. 3) are so connected that compensation occurs when the relay is energized, 





When the input signal is not changing (i.e., = 0), the circuit functions as follows. As C, charges 


If | dUjn/dt| # 0, the grid potential of L, increases, Ug decreases, L, is cut off via capacitor C,, and 
the relay is inhibited from operating until either dUjn/dt becomes close to zero or sufficient time passes for 
capacitor C, to overcharge to the new value of Ug. 


Thus, the circuit carries out a logical operation: to choose moments of time for drift compensation which 
guarantee minimum error of measurement. 


In the circuit just considered, as in the previous one, overcharging of the input capacitor, and the deleterious 
consequences connected with this, were eliminated. 


It should be mentioned that, in a number of cases, 
particularly when maximum sensitivity is required, com- 
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d; 
—— “f)| a satisf ts are provided by designs with 
ates * a() pening an + sadume vonpenaeieaa in Ri Md ond. This 

+ is explained by the circumstance that, with very high 
current sensitivity, there is an unavoidable significant 
Fig. 5. Design for an input circuit for measuring time constant at the input which limits the speed of vari~ 
small currents: a) without compensation of input ation of Ujp and thereby lowers the error deriving from 
capacitance; b) with input capacitance compen- the short-term switching off of the input signal. In fact, 
sation. suppose it is necessary to attain a sensitivity of 1 - 10~** 
amp , with a signal-to-noise ratio equal to unity, and 






















with an input impedance of not more than 10 ohms (a larger value of R would be less stable), Two cases are 
possible: a) when the input network consists of a high-ohmic resistance, shunted by an input capacitor (Fig. 5a); 
b) when the input impedance has only active components (Fig. 5b). 


As is well known [3], the signal-to-noise ratio in the circuit of Fig. 5a, using a dynamic condenser, equals 


J, ,R 
j= a ual > (2) 
kT 
C 
where k = 1,38 - 10°™ erg / degree is Boltzmann's constant, T is the absolute temperature of the resistor R, C is 
the equivalent input capacity,and Ij is the current to be measured. 


The necessary value of Tin, the input time constant, we find by substituting C = Tin/ R in (2): 


wt (3) 


ies ee XT 





By the conditions of the problem, 9 = 1, lin = 107** amp , and R = 10" ohms. 
Substituting these values in (3), and letting T = 300°, we get Tip = 4.2 seconds. 
If C is artificially compensated [8], we then obtain the circuit of Fig. 5b, for which 


Ii, R 


—S an ? (4) 
V 4kT RA} 


j= 


where Af is the amplifier's pass band which, in order to provide the given sensitivity, must have the value 
Af = —2—. (5) 


Substituting the values assumed above in (5), we obtain A f = 0.18 cycles, which is equivalent to an en- 
velope time constant of about 2 seconds. With suddenly changing signals, the time necessary to establish the 
output voltage is, in practice, about four or five times Tin, i.e.,approximately 20 seconds in the example given. 
Obviously, switching off the input signal for 0.1-0.2 seconds for drift compensation cannot generate a significant 
error of measurement. For such reasons, the circuit for the amplifier developed by us, which provides maximum 
sensitivity, was not designed to afford control of the pulse generator‘s operation. 


2. Electrometric Amplifiers with Automatic Zero Drift Compensation 





Based on the investigation discussed, two designs were developed for instruments for the measurement of 
very small, slowly varying currents, and for emf's with large internal impedances. 


We give here a short description of the type ESU-1 self-balancing electrometric amplifier. The block 
schematic of an analogous device was given in Fig. 3, but here there is provided no control over the pulse gener- 
ator's period. Fig. 6 gives the main details of the circuit, 


The instrument consists of the following assemblies: modulator-dynamic condenser C,, input stage (tube 
L), two-stage narrow-band voltage amplifier (one 6N2P tube), a power amplifier (L,), a demodulator (Tr,, 
D, — Dg, Ryo — Ryg ), an Output cathode follower (L,;), memory condensers C, and Cg, control relay P;, a carrier - 
frequency generator (to supply the dynamic condenser and the demodulator), a pulse generator (for control relay 
P,), and a supply block. Both generators are contained in one tube envelope (a 6N1P). The supply contains a 
DG-Ts27 rectifier, a filter and two stabilized voltage sources. All tubes use an ac filament supply. 


The circuit works in the following way. The input signal is divided by the removable high-ohmic resistor 
Rj. Depending on the range of current measurement, the valve of Rj may be chosen between 10° and 2 - 10" 
ohms. For measuring emf's, Ry can generally be removed, The signal voltage across R, is supplied to the plates 
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of the dynamic condenser and is modulated. The ac component across C, and C; is applied to the grids of the 
input tubes, Tube L, has a grid current on the order of 10°"! amp , which engenders a voltage drop across Rg. 
When this grid current varies, current flows across C, and C;, which can be modulated and can generate zero 
drift. However, changes of grid current in L, occur slowly, and the time constants C,R, and C3R, are sufficiently 
small so that zero drift due to the action of grid current does not exceed the noise level. 


The amplified voltage goes to the demodulator, the output of which, via the Rg-C, filter, is connected to 
the grid of the output stage L,, Part of the output voltage, picked off dividers Ryg-Ryg, is supplied as negative 
feedback through Cy to the lower plate of the dynamic condenser. 


Once a minute relay P, operates, switching contacts P,,, Pyp and Py¢ to the lower state. With this, mem- 
ory condenser C, supplies to the output the voltage existing just prior to the moment of switching, the input sig- 
nal is cut off,and the amplifier (exclusive of the output stage) is shunted, through Ry, by a 100% negative feed- 
back path (analogous to the schematic in Fig. 3). 


Condenser Cy is connected to the demodulator output and is charged up to the voltage corresponding to 
the null shift. After this, the circuit is returned to the initial state, but the null shift which had occurred is, for 
all practical purposes, compensated, 


Input relay contact Py, executes separately all the measures considered above for lowering errors and is 
mechanically connected with the armature of a standard type RKN relay. In addition to input terminals 1, 2 
and output terminals 4,5, the device has a terminal 3 on which the input voltage is repeated. Joining the shield- 
ing input lead to this terminal decreases sharply the input time constant. 


The demodulator, transformer,and carrier generator output transformer are assembled in one block with a 
special system of shielding which inhibits the appearance of parasitic voltages with respect to ground at the de- 
modulator output. This latter is very important, since during drift compensation the demodulator output is not 
directly grounded, and is connected to the amplifier input. 


A model of the instrument was constructed as one block together with its supply, and has dimensions of 
210 x 155 x 145 mm. 


Testing of the ESU-1 showed that its sensitivity is limited only by the thermal noise in the input circuit. 
Zero drift during a virtually arbitrary time interval (the observations were carried out over a two-month period) 
did not exceed twice the noise level when the internal impedance of the signal source was 10° ohms, and even 
10" ohms. For Rj = 104 ohms and the input capacitance equal to 30 micromicrofarads, the theoretical value of 
thermal noise is 12 microvolts. The experimentally determined mean-square value of the noise was 15 micro- 
volts, Thus, for all practical purposes, the theoretical limit of sensitivity was attained. The error in voltage 
measurement when the internal impedance of the source was 10° ohms was not more than +0.1% of the quantity 
to be measured + 20 microvolts. Here, we do not include the errors of the indicating or registering devices hooked 
to the output of the ESU-1. If a device of class 0.1, with a 2-volt scale, is used, the ESU-1 provides the follow- 
ing limits on measurement: 


1) 0 = 2 millivolts with an error of t 1% + 20 microvolts; 
2) 0 - 20 millivolts with an error of + 0.2%; 
3) 0 = 200 millivolts with an error of + 0.2%. 


The current sensitivity is proportional to the size of the resistance Rj, through which the current to be mea- 
sured is passed, For Rj = 2 * 10’? ohms, the current sensitivity equals 1.- 107" amp . With this, the zero drift is 
determined, not by variations of the contact emf inside the dynamic condenser, but by the currents in the input 
insulators, and by ionizations currents created by cosmic radiation and by radioactive contamination.The magnitude 
of the current drift can be lowered to 5 * 107"? amp. during an 8-hour period. A decrease in the size of Rj pro- 
portionately decreases the current sensitivity and the input time constant (Tip = RjCjn). For Rj = 4° 10° ohms, 
the current sensitivity is 5 + 107" amp., and Tj, = 0.1 seconds. The error in current measurement is almost com- 
pletely determined by the instability of the high-ohmic resistor Rj, where the error increases as Rj increases, 
Thanks to the high sensitivity of the ESU-1, and the stability of its null, the required value of Rj is much lower 
than is needed for measuring the same currents by other instruments, which significantly increases the accuracy 
in measuring small currents. Table 2 gives the basic parameters of the ESU-1, and also those of the best domestic 
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TABLE 2 











Electrometer type 








EMU-2 (HII 

Parameters and Characteristics ESU-1 broadcasting “Cary model | "Vibron elec- 
receiver and 31°° trometer model 
acoustic, 33b"* « 
MRTP SSSR) 

Number of envelopes 4 1 —*es 4 


Dimensions in millimeters 


Voltage sensitivity in millivolts 


Current sensitivity in amperes, 
for Ry = 2° 10” ohms. 


Input impedance in ohms 


Input capacitance in micro- 
microfarads 


Daily zero drift in millivolts 
Output impedance in ohms 


Maximum output current in 
milliamperes 


Supply 





210 x 155 x 145 


0.02 


1x 107!" 
10%* 


30 
0.02 
800 


1 


220 volt network 
without a stabili- 
zer 





*Developed by “Applied Physics Corporation” (USA). 
* *Developed by “Electronic Instruments Ltd." (England). 


* **Not provided in the sketches, 


260 x 225 x 215 


(without the supply) 


10 


2x 10° 


10” 


VS-14 stabilizer 





0.02 


1x 107" 


25 
#0.2 


without 
stabilizer 


380 x 290 x 200 


0.1 


5x 107% 
10" 


£0.2 
1500 


1 


110 volt network, 
220 volts with- 
out a stabilizer 





and foreign electrometers. It is clear from the table that our instrument excels, in a number of parameters, 


similar 


electrometers of other types, in particular, in its daily zero drift, which is 10 to 50 times lower than 


that allowed by the other instruments. When drift is observed over even longer periods of time, the ESU-1 shows 


itself to be even more superior, 


SUMMARY 


1. As a result of an investigation into the possibility of zero drift compensation in electrometric ampli- 
fiers, the sources of errors arising in compensation were explained, and methods of lowering them were suggested. 


2. Tests of a model of the instrument developed attested to the superiority of employing zero drift com~- 
pensation, in amplifiers with dynamic condensers, by means of contact switching. 
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ELECTRIC METHODS OF FREQUENCY CONTROL FOR STABLE RC-GENERATORS 


V.P. Demeshin 


(Moscow) 


A new electric method is considered which, within wide limits, allows the 
frequency of stable RC-generators to be controlled, Results are presented of in- 
vestigations into the dependence of frequency on the controlling voltage and into 
the character of the transient response, with both the noninertial and the inertial 
nonlinearities being taken into account, Practical circuits for controlling RC- 
generators are given, 


Electric methods of controlling oscillation frequency in generators, although widely used in communica~ 
tion work do not, as a rule, satisfy the requirements on stability and deviations imposed in telemetry. Today 
there have been developed comparatively stable frequency-type transformers which are the bases of generators 
of nonsinusoidal oscillations which are to be controlled by electrical means, said transformers functioning at 
audio and subaudio frequencies. 






































y 
f 0 a - t 
Y¥,sinwt = Ri, ' i 
Ri} U, 
2 I} r 
j a —_4 
- = iG ~ t 
-—Y, ——+ z-6 
Fig. 1. Circuit of the controlling network. Fig. 2 


In the present work we consider, following the proposal of V. A. Il"in, a simple electrical means of con- 
trolling the frequency of a high-stability RC-generator, which allows one to obtain a large variation in frequency, 
ranging from the tens of cycles to several magacycles, 


The control of frequency of an RC-generator is based on the use of circuits with elements possessing piece- 
wise linear characteristics, which are acted upon by both ac and dc voltages, 


The controlling device is a simple network, consisting of two resistors and a diode (Fig, 1). An ac volt- 
age U;psinwt is applied to the input of the network, and a dc voltage U, to the diode, The diagram of the 
voltage across the diode is given in Fig. 2. Under the condition Um > Up, the diode periodically conducts 





and cuts off. The time r during which the diode conducts during each period Tg is determined by the follow- 
ing relationship between the amplitude Uy of the ac voltage and the magnitude U, of the de voltage: 


re sin Us 
a 
Leet (1) 


Tr 





As is obvious from (1), by varying the magnitude of the dc (controlling) voltage it is possible to obtain 
various relationships between the durations of the conducting and cut-off states of the diode. 
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Fig. 3 


The conducting diode shunts resistor R,, due to which the input impedance on the side of the ac voltage 
is lowered. If one considers the average value of the network impedance over one period, then its magnitude 
will depend on the time during which the diode conducts and, consequently, on the controlling voltage. This 
dependence of the input impedance on the voltage will also be used in the. control of the frequency of an RC- 
generator, the circuit of which is given in Fig. 3. Here, to the diode is automatically applied an ac voltage, 
Umsin wt of the frequency to be generated. Conducting under the effect of the ac voltage, the diode periodic- 
ally shunts resistor Rg, as a result of which the oscillation frequency of the generator varies. If a controlling volt- 
age is applied to the diode and its magnitude is varied, it becomes possible to control the frequency of the oscil- 
lations to be generated. 


The actual frequency of oscillation in a generator with a diode controlling device can be given as the 
average of its values during the period 


@ Tt + @, (T — T) 


o> T 





By wg and w, here should be understood the equivalent frequencies of the oscillatory processes for the times 
when, respectively, the diode conducts and the diode is cut off. 


Oscillatory Processes in the Simplest Transformer with One Diode 


The dependence of the frequency of oscillation on the controlling voltage, and other characteristics of the 
transformer, are determined, to a large extent, by the type of generator. We shall limit ourselves to the consider- 
ation of one type of high-stability controlled RC-generator with a parallel-series phase-shifting circuit (with a 
Wien bridge) [1]. Its block schematic and its equivalent circuit are shown, respectively, in Fig. 4 and Fig. 5. 
The controlling diode can be connected in the phase-shifting circuit of the generator in various ways. The simp- 
lest of these ‘is shown in Fig, 6. We shall consider in detail the properties of a generator with one controlling 
diode in the Rg circuit in order to establish the fundamental regularities. 
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The electrical processes in a generator with a parallel~series 




















+5 phase-shifting circuit are described by a second-order differential equa- 
R, | l tion [2] 
‘ t U + 280 + wv =0. (2) 
U, 
2 F |, he The solution of this equation has an oscillatory character, and 
4 mee" 4% may be presented in the form 
) U (t) = Ume— sin (wt + 9), 
~ j 
J 
“ .'? (ea. K 
Fig. 6 8 = > eh (et ta + 3 ee 5K) ; (3) 


of aii anes 
RikC,C, * 

The introduction of a controlling diode in the generator circuit has as a consequence that the resistance 
R, in the phase-shifting circuit periodically changes its value. In accordance with (3), this has the result that 
the frequency of oscillation and the damping coefficient also vary. The influence of the diode significantly 
complicates the process in the generator, since the system amplifier~phase~-shifting circuit with diode acquires 
a sharply expressed nonlinear character. In addition to the noninertial nonlinearity, the source of which is the 
diode, an inertial nonlinearity is introduced into the generator design in the form of a thermal resistance (ther~ 
mistor) in the negative feedback path [2]. Each of these nonlinearities has its effect on the oscillatory process. 
The noninertial nonlinearity shows its influence principally on the form of the voltage to be generated and on 
its period of oscillation; the inertial nonlinearity determines the amplitude of the oscillation and the character 
of its establishment. 


We now determine the frequency of oscillation in a generator with a controlling diode, and the dependence 
of this frequency on the controlling voltage. To do this, we first consider the form of the oscillation in the steady 
state without considering the influence of the inertial nonlinearity, since its time constant is generally chosen to 
be many times greater than the period of oscillation. 


The differential equation of an RC-generator with a parallel-series phase-shifting circuit and a controlling 
diode may be written in the following form: 


U +n + w =0, (4) 
where 
_<> 5) fue filgoy Gls K i 
n= Ror (or tam +1— T= pe) = Tce 


The coefficient a is a function of the ratio of the controlling voltage to the instantaneous value of the ac 
voltage on the diode. Numerically it equals one when Uppsinwt/U, <1, and equals the ratio R$ Ay when 








Umsinet/U,>1. Here R§ and Rj are the resistances in the R, network when, respectively, the diode is con- 
ducting and is cut off. The change in the coefficient a occurs abruptly at the moment when Umsinut = Uy. 
According to (4), an abrupt change of a entails changes in the frequency of oscillation and in the damping coef- 
ficient 6 = n/2. Thus, the coefficients of Eq. (4) are functions of an independent variable, and are changed 
abruptly at a definite value of this variable. 


The nature of the voltage variation across condenser C, can be determined by successive integrations of 
Eq. (4). For this, we divide the period of oscillation into two parts, one of these corresponding to the time when 
the diode is cut off, the other to the time when the diode is conducting. In each of these segments Eq. (4) will 
be considered as linear, and with constant coefficients. At the transition points the solutions will be “spliced” 
by taking into account the initial conditions. 


We denote by n, and n, the corresponding values of the coefficient n for a cut-off and a conducting diode, 
Then 


a oe. K 
m= oat at! — ier): 

ti. 1 C; Ri cS ae 5 
"= aor (ot tam t+!— tar): ” 


For RC-generators with parallel-series phase-shifting circuits, one generally chooses Ry = Rg = R, Cy = Cy = 
=C,K >>1[3]. The ratio R$/R} we take equal to 0.5. This gives the optimal relationship between the fre- 
quency deviation and the attendant amplitude modulation. With these relationships taken into account, we get 


m = Re (3+): m = Ro (4— >): (6) 


Assuming that the negative feedback coefficient 8 remains constant for a small number of periods, we ob- 
tain from (6) 


or. ae. (6a) 


Without going through the intermediate steps of the calculation, we write just the final expression for the 
solution of the equation for an RC-generator with a parallel-series phase-shifting circuit; 


U (t) = ((U (0) n + 20 (0)| sin at + 20 (0) cos wt} + e—*, (1) 


where 





x n 4 4 n 
=F. om ge Vi-F- 
By introducing the notation b, = [U (0)n + 2U (0)) we and by = U(0), we can write Eq. (7) in a simpler 
form: 
U (t) = Um, Be sin (wt + 9), (8) 


where 


B=VU+B, o=arctan — 


The solution (8) was obtained for an arbitrary value of the coefficient n and is therefore valid both for the 
conducting and for the cut-off diode. 


We now follow through the course of the oscillatory process, starting at the moment when the diode begins 
to conduct. For this, we introduce the initial couditions into the solution (8). As one of the initial conditions 
might be taken the voltage on condenser C, 








U (0) = Une sin wt, (9) 


where w, and 6, are, respectively, the frequency and damping coefficient with the diode cut off, and t, is the 
time when the diode begins conducting. 


As the second initial condition, we choose the value of the derivative of the voltage on condenser C,. To 
determine the derivative, we consider the current distribution in the phase-shifting circuit with the diode (Fig.6). 
We set up the equations for the current at node A for two moments of time: immediately before the diode be- 
gins to conduct, and immediately after: 


L=h+h= Ua + Co, (10) 
L=h+ hy =Uizge + Cl (11) 


Here 1}, U4 and Ij, Uj are the currents and voltages , respectively, up to and after the moment at which 
the diode begins to conduct. The voltage on condenser C, cannot change instantaneously at the moment of 
diode switching and, therefore, neither can the current I,. With this in mind, we subtract both members of (11) 
from the corresponding membezs of (10). As a result, we obtain 


Ui + Cll, — Ue — Ci we (12) 


or 
** a i /i—a 
0, = 01 — a0 (SS) 
Setting a = 0.5, we obtain the value of the derivative of the voltage on the condenser at the moment when 
the condenser begins to conduct: 





U* (0) = 0, —eU;. (13) 


We now introduce initial conditions (9) and (13) into Expression (8) for the voltage on the condenser. After 
grouping terms,we obtain 


Uy (t) = U mes {{(s 4+ 28, — 2 a) sin «,t + 2e0 608 wt, | e— sin ant + 


(14) 
+ 20,e—* sin at; £08 «yt } —— : 
We replace 6 by n,/2 and also, from (6a), we express n, in terms of ng, obtaining 
U (t) = Ue {{[ 2eo, 008 cyt, — sin coyt,) sin cogt + ane 


++ 20d, Sin wt; COS Wat} ae. 
Denoting the amplitude of the oscillation Um,e~®1% by Um, and also introducing the additional notation 


b= te (2; Cos o,f, — Sin@t,), by = sina,t,, (16) 


we get 
U,(t) = BLU me sin (ast + Pa), (17) 


where 


By=Vii+ti, ge =arctan Z. 
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| i m | i e 
Ng=0.2 @y=1.41 RC ; g=0.4 wii? Ro ; %=0.6 w=! 38 Re : m_—0.8 a ; 
t 
°° Bal = 0.44 St = 0.93 5- St = 1.36 St = 1,86 —> 
B, o B, o By m By o 
0 0.65 0 0.68 0 0.71 0 0.738 0 
10 0.6 16°30’ 0.632 15° 0.665 15° 0.68 15° 
20 0.598 35° 0.62 33°50’ 0.643 32° 0.66 31° 
30 0.63 52°30’ 0.646 50°30’ 0.66 48° 0.675 48° 
40 0.692 68° 0.705 65°30’ 0.726 62° 0.72 63° 
50 0.78 79° 0.79 78° 0.8 77° 0.78 77° 
60 0,867 90° 0.865 90° 0.87 90° 0.865 90° 
70 0.94 97° 0.94 97° 0.945 96° 0.94 96° 
80 1.02 104° 1.02 104° 1.02 104° i 103° 
90 1.05 110° 1.06 110° 1.06 110° 1.065 109° 


As is clear from (17), after the diode begins conducting ,the process retains its oscillatory character, but 
with different amplitude, phase,and damping coefficient. The new values of these parameters will depend on 
Ng and on the initial phase gy, = w,t, with which the diode begins conducting. The initial phase g, is deter- 
mined by the magnitude of the controlling voltage, but the damping coefficient is established under the action 
of the inertial negative feedback in such a way that the steady-state conditions are met. There has been no 
success in finding analytic conditions for the essentials of the steady-state oscillations. Therefore, they are de- 
termined by the method of successive approximations. In Table 1 are given the parameters of the process after 
cutting off of the diode, as they depend on the initial phase (magnitude of the controlling voltage), for succes- 
sive series of values of the coefficient n,. 


For each value of the coefficient ng a complete cycle of oscillation was defined and its steady-state 
character verified. 


The diode is cut off again after the instantaneous total voltage on it becomes less than zero, The character 
of the voltage variation across condenser C, with the diode cut off is determined by Eq. (4) with the initial con- 
ditions being taken into account. For the initial conditions, as previously, we take the voltage on condenser C, 


U,(0) = Umsin ¢ (18) 


and its derivative, To determine the derivative, we consider the current distribution at node A (Fig. 7). With 
this, it should be borne in mind that the current I, is shifted in phase by less than 4/2 with respect to the volt- 
age on the condenser (Fig. 8). Due to this, I, can have either a positive or a negative voltage, depending on 

the initial phase. If current flows to the node then, due to the decrease of I,, the derivative of the voltage across 
the condenser after cut-off of the diode may change sign. However, in practice the voltage on the condenser 
cannot increase due to the limiting action of the diode, As a result, something akin to a sliding regimen will 
ensue. Since the diode, as a relay element, has a dead zone close to zero, the sliding time is also close to zero. 
The voltage on the condenser only begins to decrease after the following condition is fulfilled: 


ys Ih. (19) 
We now write the equation for the current at the node after the diode is cut off; 
1, = Ig — Is = 1, — CUg. (20) 


Introducing into (20) Condition (19), we get the initial value of the derivative of the voltage on the con- 
denser at the moment of diode cut-off (cf. Table): 
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TABLE 2 U"(0) = 0, (21) 
With an initial phase such that ,, at the moment 
ot at di boagel Pe af % of cut-off, has a negative value, the derivative of the 


voltage on the condenser changes its value without 
changing sign. The equations for the currents at node 
A, in this case, are written in the following form: 




















0 0.53 | 4.30] 180° | 44.0 

0.174 | 40° | 0.60 | 4 157° | 42.0 

0.342 | 20° | 0.671 0.76 | 102°| 9:3 —1I,—I,+/,=0, 

0.500 | 30° | 0.74 10.97] 97° 6.9 I r ° (22) 
0.642 | 40° | 0.80 | 4 96° | 5.0 —4,—1,+1,;=0. 

0. 66 | 50° |0.88/14.40] 95°} 4.5 


Subtracting, and expressing the currents in terms 
of voltage, we get 


ae ee * U. © 
Ui, = 0, 0)+—te (23) 


Taking initial conditions (18), (21) and (23) into account, we obtain from (8) the corresponding two ex- 
pressions for the voltage on the condenser: 


U;,(t) = By By mee sin (at + Gn), Bi = V bi, + bis, 


1 P A bis (24) 
by, = om SiN Wate, big = SiN @st, Ge, =arctan “he 
U} (t) = BiB U mere sin (yt + G12), By = V 0%, + bbe, 
bes = SiN (Wat + 92), 5g, = — [Sin (cogt, + 9g) + Qed, COS (Waly + 92)], (25) 
Pie =arc tan e . 


By employing formulas (24) and (25), one can determine the complete oscillatory cycle in a generator 
with one controlling diode. 





$f, % In Table 2 are given the computed data which 
: characterize the process, and the coefficient n, in the 
t steady state, 


* Fig. 9 shows the dependence of the relative fre- 


quency variation on the ratio of the controlling voltage 
to the amplitude of the ac voltage on the diode, Curve 
2 was obtained by computation, and curve 1 was ob- 
tained experimentally. The insignificant amount of 
the discrepancy between them is evidence of the cor- 
rectness of the results we have obtained. 


In a bridge-type generator with a parallel-series 


se.0 phase-shifting circuit and one controlling diode, the 
8. maximum frequency deviation reaches 14%, However, 








the portion on which the dependence of the frequency on the controlling voltage deviates from linearity by 
more than 4% does not exceed 8%, Roughly, with an error not greater than 20%, this dependence can be com- 
puted from the formula 


_ Us 
are 3in 7; 
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Establishing the Steady State in a Generator with a Controlling Diode 





The character of the transient response is very important in evaluating a transducer, since this determines 
a very important characteristic of the transducer — its speed of action. 


The analysis of the transient response in the generator we have been considering can be reduced to solving 
a system of two equations, the first of which describes the generator without considering inertia, and the second 
of which describes the properties of the negative feedback loop with nonlinear inertial elements in it: 


0+4+U=F(U,U,0, 8, p, 6), B =7(0, B), (27) 


where U is the voltage on the first tube's grid and 6 is the coefficient of the negative feedback path. 
In setting up Eq. (27), the influence of the controlling diode was factored in by approximating the char- 
acteristics of the amplifier and phase-shifting circuits by a fifth-degree polynomial 
I = SU(U + jU* + yu® + zu. (28) 
The thermistor characteristics can be presented as a power series of the following form: 


p = py (1 + a8 + b0? + 20%), (29) 


where p, is the filament impedance at zero temperature. 
For tungsten the coefficients in Formula (29) have the following values; a = 4.6 - 107° degree ~*; b = 
= 0.7°10°' degree “*; g = 0,062 °10°* degree ~°. 


We obtain an approximate solution to System (27) by using the well-known method of small parameters 
for the oscillation amplitude, in the following form: 


O=*(—F4(14+ Z)+(1—Gyr+ gv) — 38H], (30) 


where k is the amplifier gain, 8 is the coefficient of the negative feedback path, y and z are coefficients of the 
series which approximates the amplifier characteristics with the diode being taken into account. 


By setting UW equal to zero in (30), and remembering that k >> 1, we obtain the equations which determine 
the steady-state value of the oscillation amplitude: 


3 5 
U (1— fA? + g2A*) =0, B=0. (31) 


The fact that there are two equations for the steady-state value of the amplitude is explained by the pres- 
ence in the circuit of the investigated generator of two nonlinear elements, noninertial and inertial. The ampli- 
tude of oscillation after the steady state is established can be presented as the sum Ugg = Ugy + AU, where Ugy 
is the value of the oscillation amplitude obtained from (31) and Ug is the steady-state value of the oscillation 
amplitude. The size of the correction is found by solving, simultaneously with Eq. (31), the equation 









































Fig. 11 
AU = U,, (1 — e-® sin (Q# + 9)], 
br =a [(1—2) + (Fy As +-g-243) AC — By], (32) 


where T is the thermistor time constant and A is a coefficient depending on the amplitude of the oscillation. 


As is obvious from (32), the process of establishing the steady-state amplitude has a damped oscillatory 
character, Computations, carried out for an actual generator system, show that the damping is basically deter- 


mined by the thermistor time constant 6 7 = > (0.5 + 0.4— 0,12) = 0.78 > - Depending on the thermistor 


type, the time constant can vary within broad limits — from tenths of a second to several seconds, Fig. 10 gives 
the oscillograms for the process of establishing the steady-state amplitude for generators with various types of 
thermistors. The oscillograms of Fig. 10a and b were obtained for a type TP*/,g thermistor, where there was a 
larger negative feedback in the case of oscillogram a. When the negative feedback was decreased, the damping 








was almost not decreased at all, which is in complete correspondence with the results we have obtained. In a 
generator with a type TP*/ os thermistor, the time needed to establish the steady state is approximately 0.5 sec- 
ond . The oscillogram of Fig. 10c was obtained with a type M-28 incandescent tube in the role of the ther- 
mistor. The oscillogram of Fig. 10d was obtained with a type T-9 beaded semiconductor thermistor. The time 
taken to establish the steady state was about 0.1 second. 


In using thermistors in generator designs, it should be kept in mind that the working temperatures of the 
majority of semiconductor thermistors lie within room temperature limits [5], and variation of the ambient 
temperature leads to a variation of the amplitude of oscillation. Therefore, for transformers of the highest 
stability, thermistors with high-temperature working points should be chosen. 


The results of analysis of the simplest circuits for controlling the frequency of RC-generators with a diode 
in one portion of a phase-shifting circuit allow the following conclusions to be drawn; 


1, Frequency control of an RC-generator can be implemented over the entire frequency range in which 
such generators are capable of functioning,i.e., in the range from several cycles to several magacycles. 


2. If the elements of the phase-shifting circuit are correctly chosen, instability in the parameters of the 
controlling diode has only a small influence on the generator frequency. 


3. A transformer with one controlling diode allows frequency variation within approximately 14% limits 
to be obtained with an attendant amplitude modulation of 40%. The amplitude modulation has as a consequence 
that the damping coefficient depends on the controlling voltage. 


4. A decrease in the accompanying amplitude modulation and an increase in the frequency deviation can 
be attained by introducing controlling diodes in both sections of the phase-shifting circuit. 


5. The speed of action of the transformer is basically determined by the thermistor time constant. 


Transformer Circuits with Controlling Diodes 





The results obtained from the investigation of a generator with one controlling diode allowed a number of 
other transformer systems to be developed, systems in which the same control principles were used, but in which 
a number of important shortcomings, inherent in the simplest systems, were avoided. 


Fig. 11 shows the circuit of an RC-generator in which diodes were connected in both portions of the phase~- 
shifting network, With such connections, the magnitude of the resistance is lowered in both portions of the bridge, 
as a result of which the effective oscillation frequency deviation is increased and the accompanying amplitude 
modulation is decreased. These results correspond completely to the conclusions previously drawn, According 
to (4), if diodes are introduced into both portions of the phase-shifting network then, in the expression for the 
oscillation frequency, the coefficient a is squared, and in the expression for the damping coefficient, the coeffi- 
cient a is cancelled out, This is valid : only in the case when the currents through the impedance R,R, coincide 
in phase, Actually, the currents are shifted in phase, and the diodes do not begin to conduct simultaneously, 

i.e., there is, with this, an interval of time when only one of the diodes conducts, which leads to a decrease of 
the frequency deviation and an increase in the accompanying amplitude modulation. 


Fig. 12 shows the dependence of the frequency and amplitude of oscillation on the controlling voltage, 
with the ratios R[/R} = R%/R} = 0.3. The fundamental disadvantages of this circuit are the significant amount 
of amplitude modulation, and the presence of additional networks,impairing the bridge properties. 


Fig. 13 shows another circuit for controlling an RC-generator with one diode. In this circuit, the diode is 
connected between two phase-shifting networks of the same type but with different time constants. The tube's 
grid circuit has a large time constant. The circuit elements are so chosen that the transfer constant does not 
change when the diode resistance drops, Under the action of the ac and controlling voltages, the diode periodi 
cally conducts and shuts off, and first one, then both portions of the phase-shifting circuit are periodically switched 
into the grid circuit. The oscillation frequency of the generator varies in dependence on the ratio of the times 
in the conducting and cut-off states of the diode. By choosing network portions with different time constants, 
one can, within broad limits, vary the generator's oscillation frequency, 


Fig. 14 shows the dependence of the oscillation frequency on the controlling voltage when the ratio of 
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af, network time constants equals six. The oscillation frequency changes 
7* by more than a factor of three. The portion with a nonlinearity less 
aoa} than 2% attains 75% of the entire range of frequency variation. Fig. 
15 gives curves characterizing the magnitude of the maximum frequency 
a+ deviation with control for different values of the network time constant 
ie A rae ratio. If the ratio of time constants is greater than ten, the frequency 
“20-0 0 2 re varies by a factor of six and more, Transformers with such a method of 
1 controlling the frequency of an RC-generator with a parallel-series 
phase-shifting network possess high stability. Fig. 16 shows the depend- 
ence of frequency on the voltage on the tubes anode, If the supply volt- 
age varies from —15% to +35%, the frequency varies only by 0.08%, 
Fig. 16 and if the supply voltage varies by 1%, the frequency varies only by 
0.004%. An inherent, unavoidable disadvantage of a transformer with 
a controlling diode is the distortion in the form of the generated voltage. In Fig. 14 is given the curve of the 


magnitude of the second harmonic (the largest harmonic) in the generated voltage as a function of the control- 
ling voltage. 

















In conclusion, it should be mentioned that the electrical method considered here allows the harmonic 
oscillatory frequency of a high-stability RC-generator to be controlled within limits not attainable by other 
known methods of control, Its stability satisfies the requirements imposed by telemetry technology, and it may 
be used in a very wide frequency range, since the limitations are due only to the generator properties. The trans- 
former with diode control is simple, reliable, and does not require complex tuning. Similar methods of frequency 
control can find use in telemetry and also in measuring technology and in other branches of technology. 
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ANALYSIS OF VARIOUS DESIGNS FOR BROAD-BAND TUNING 
OF HIGH-FREQUENCY WAVE-TRAPS 


G. V. Mikutskii 


(Moscow ) 


Various schemes are considered for broad-band tuning of high-frequency wave- 
traps. Generalized curves are given for facilitating engineering computations of the 


tuning elements, and recommendations are given regarding the protection of these 
elements from overvoltage. 


In high-frequency electrical transmission channels using wires, such as telephone communications, remote 
control and relays, high-frequency wavertraps are used for protective purposes, namely, to inhibit the flow of 
high-frequency currents along the substation busses. The basic element of a wavestrap is a power coil; the 
higher the inductance of the power coil, the better the trapping properties of the wave-trap. However, the trap- 
ping powers are determined, not only by the inductance of the power coil, but also by the design for tuning this 
coil. An analysis is provided below of various designs for broad-band wave-trap tuning, and a comparison Is 
made of broad-band tuning with ordinary single-frequency tuning. 


1. Loss Introduced by the Wave-Trap 





The loss introduced by a wave-trap is defined, for the circuit of Fig. 1, as one-half the natural logarithm 
of the ratio of the power produced by the line generator for Zy = @, P;, to the power P,, put out by the generator 
with a finite value of impedance Zy (| Zw| #0) [1]. Starting from 


this, the loss introduced by the wave-trap is determined from the expres- 
sion 





= Inji+ 


lo 


(1) 





Z in 
——7 |» 
2,(1+32) 
Fig. 1. Circuit for computing the 


loss introduced by a wave-trap. where Zip is the input impedance of the electrical transmission line, 

Zw is the series-connected wave-trap impedance and the substdtion's 
internal impedance, and Zg is the internal impedance of the equivalent generator connected in the line at the 
point where the capacitor coupling is effected. 


If it is assumed that the generator is matched with the line, which has purely active input impedance, then 
Zg = Zin = tn- With this, 


bo = In| 1+ 38 ‘ (2) 





In the general case, impedance Zy is complex, 1.e., Zw = tw + iXw. If Zy is purely active, then Eq. (2) 
takes the form 
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b= In(1 + zn), (3) 


2r 
Ww 
and, for small imaginary Zw (purely reactive), 
b, me tin( 44 "in (3a) 
lo 2 4x? | * 
w 


Analysis of Eqs. (2), (3) and (3a) allows the following conclusions to be drawn; 


8) Active impedance has a significantly greater bypassing action on the channel than does reactive. For 
example, for tw = tn (Xw = 0), the insertion loss is 0.41 nepers while for xw = tp (tw = 0), the insertion loss is 
0.11 nepers. Therefore, if only the reactive component of the total wave~-trap impedance is used, this latter can 
be significantly smaller than the magnitude of the active impedance if only the active impedance were to be 
used to obtain the same amount of insertion loss. For an insertion loss of 0.41 nepers, it would suffice to have the 
reactive impedance equal only to 0.45 rp. 


In the overwhelming majority of the cases, the input impedance of the substation has a capacitive character, 
Therefore, there is a real danger of compensating the reactive impedance of the wave-trap if the latter has an 
inductive character. If the wave-trap impedance is capacitive in nature, then such compensation is unlikely, 
and can be neglected in the majority of cases for telephone and remote control applications. 


b) If the insertion loss is determined by the reactive impedance, then adding an active impedance in series 
with it somewhat increases the insertion loss. For example, for ry = 0 and Xw = %p, bjo = 0.11 nepers, while for 


tw = Mn and Xw = tin, big = 0.24 nepers. 


c) If the insertion loss is only determined by active impedance, then adding a reactive impedance in series 
with it, no matter what the sign of the impedance, will lower the insertion loss, Therefore, for high-frequency 
relay protection channels, as the most susceptible, the orientation should be towards only the active component 
of the total wave-trap impedance. 


The input impedance of the substation, connected in series with the wave-trap, may in some cases lower 
essentially the magnitude of the loss introduced by the wave-trap. However, the magnitude of the input imped- 
ance of the substation is not constant, It depends on the amount of equipment connected to the substation, the 
switching of this equipment, etc. Therefore, in computing the loss introduced by the wave-trap, the worst case 
is generally taken, i.e., the input impedance of the substation is considered either as zero or as some magnitude 
which can completely compensate the inductive component of the total wave-trap impedance. 


2. Single-Frequency Resonant Wave-Trap 





A single-frequency wave-trap is a simple parallel resonant loop formed by the power coil and the tuning 
condenser, With a ‘sufficiently high Q-factor of the circuit, Q, = wL/r, the actual resonant circuit (Fig. 2a) can 
be replaced by a circuit without losses, shunted by a resistance r, equal to the resonant impedance of the actual 
Circuit Zres (Fig. 2b): 


tg = Zres = wmLQ, 
where wy = i is the frequency of the tuned circuit. 


The dependence of the total impedance of the circuit of Fig. 2b on the frequency can be given by the fol- 
lowing expression; 
Z WwW i u 


ra ipua—lipe (4) 





where u is the relative frequency detuning. 
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1 u=@(=—-2). (5a) 
m, tials 
a . where Qb . 
Fig. 2. Circuits of single-frequency tuned wave~-traps. The dependence on the frequency of the active 


and reactive components of the total impedance of 
the circuit of Fig. 2b is given in Fig. 3. 


For widening the frequency band within whose 
limits the active component of the total wave-trap 
impedance is sufficiently great, the resonant imped- 
ance is frequently decreased by connecting an active 
impedance rc (Fig. 2c) in the series with the tuning 
condenser. Such a wave-trap falls in the category of 
broad-band wave~traps, and received the name of 
“dulled” resonant circuit, 


The dependence of the total impedance on fre- 
quency for a dulled circuit (Fig. 2c) is determined by 


4, 





the expression 
Fig. 3. Dependence on frequency of the active and Zy 1 os 
reactive components of the total wave-trap impedance hb te a( =) ae 


in the circuit of Fig. 2b. 


s u 
ipa (: oe adie 
Computation of the dependence of the active and reactive components of the total impedance of a dulled 
circuit can be carried out using the graphs of Fig. 3 by multiplying the quantity ry/ 1, obtained from the graph, 


ay-1) 


by (w/wm)*, and the quantity xy /r, by (, = ue 


(af (8) 


The impedance r,, which enters into the expression for the relative detuning u, is determined from the ex- 
pression 


wo? L* 
o,L ° 


r. ao om (7) 


where Q is the natural Q-factor for the wave-trap circuit (for r¢ = 0), 


rs = 


3. Broad-Band Wave-Trap Based on a Band~Pass Filter Circuit 








Circuits for a broad-band wave-trap in the form of a band-pass filter are shown in Fig. 4a and b. 
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Fig. 4. Circuits for broad-band wave-traps in the form of band~pass filters, 
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The dependence of the total impedance on the frequency for the two-loop circuit of Fig. 4a is expressed 
by the formula [2] 








ou 4 .  wi—k+ ku) (8) 
mm i+0—m® e+ ee / Times ee’ 


where u is the relative frequency detuning, determined in accordance with (5), rz is the impedance of the filter 
load, k = L/13Cq, and tm = Vey, = — - 


VIC, Vic 
metric mean of the bounding frequencies of the trapping band. 





is the tuned frequency of the circuit, equal to the geo- 





For k = 0, the circuit of Fig. 4a reduces to the parallel loop circuit of Fig. 2b. The wave-trap circuit 
(Fig. 4a for k = 1) is known in the Soviet Union under the name of the Lifshits circuit [3]. 
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Fig. 5. Dependence on frequency of the active and Fig. 6. Dependence on frequency of the active and re- 
reactive terms in the total wave-trap impedance active components of the total wave-trap impedance 
using the circuit of Fig. 4a; 1) k= 0; 2) k = 0.6; using the circuit of Fig, 4b: 1) kg = 1.05 and ky = 0.55; 
3) k =1.0. 2) kg = 0.9 and kg = 0.7; 3) ky = 0.75 and kg = 0.7. 


Fig. 5 shows the dependence on frequency of the active and reactive components of the total impedance 
for the circuit of Fig.4a. The curves are symmetric with respect to the mean frequency (u = 0). Therefore, 
Fig. 5 shows only the halves of the curves which correspond to positive values of u. 


From Expression (8) can be obtained the following formulas for computing the elements of the broad-band 
wave-trap using Fig. 4a: 
{ 


L 
¢ =o ’ LL, =S e 9) 
: kr} w? Cs ( 
For the wave-trap implemented by the three-loop circuit of Fig. 4b, the dependence of the total imped- 


ance on frequency is determined from the following expression [2]: 


Zw 4 .  u(A + Bu? + Fut) 
7 1+Ce8+ Duy ew  / TC + Dut F Eu® 
Here, u is determined from Expression (5), A = 1 + ky — kg, B = ky [ky — ks (2 + kg)], C= 1— 2kg + kg(2 +s), 
D = ke [1 — 2kg(1 + k,)], E = iG , 
te gt oA 
— Vic, Vis, Vic,” 








(10) 








kg = 

















Fig. 6 gives the dependence of the active and reactive components of the total impedance on the fre- 
quency, using the circuit of Fig. 4b for various values of the coefficients ky and kg. 


Expression (10) provides the following formulas for computing the elements of a wave-trap constructed 
with the three-loop circuit: 


L 4 
we pons + L, = wc,’ 








L i 
l= +) pte (11) 


For computing the wave-trap, it is necessary that the allowable magnitude of insertion loss be given for 
the mean frequency and for the boundaries of the trapping region (band), From Formulas (2) or (3) are deter- 
mined the minimum admissible values of the active impedance at the mean frequency and the active or reac- 
tive component at the boundaries of the trapping band, From the magnitude of tw/t,_ or Xw/t,, taken from the 
graphs of Fig. 5 or Fig. 6, the relative detuning corresponding to the boundaries of the trapping band is deter- 
mined, 


If the magnitude of u is found, the boundary frequencies of the trapping band can be determined from 
Eq. (5) in the form 











ha= V fe + — fn) a fon (12) 


On the basis of (12) one determines the trapping band and the arithmetic average of the boundary fre- 
quencies of this band: 








Mf = f— hy = fy = hy 





r ry (13) 
The boundary frequencies of the trapping band will be 
h=f—, hr=fo i (15) 


By using (13) and (14) one can express the trapping band in terms of the value of the arithmetic average 
of the frequencies 





r / 
Af = a(V 1 + 4ft aey —1) as) 
With a relatively small trapping band (A f < 0.3f,), f5* fm, and Expression (16) reduces to (13), i.c., 


Af = 2nu — fp. (17) 
2 


If either the upper or lower limiting frequency of the trapping band is given, the second bounding frequency 
can be found from one of the expressions 








h= mano ae | Oe _ A= — == . (18) 
2 2 


4. Broad-Band Wave-Trap Based on a High~Pass Filter Circuit 


The circuit for a broad-band wave-trap in the form of a high-pass filter does not differ from the circuit 
for a dulled loop shown in Fig. 2b, The broad-band wave~-trap differs from the dulled loop in that the resistance 
Ic = tg is chosen to be of such magnitude that for ry =r, the loss introduced by the wave-trap does not exceed 
a permissible magnitude, The dependence of the total impedance of the wave-trap on frequency, determined by 
Expression (6), can be put into the following form: 











L ( 4 
Zw wi wl? +i oLr, — Cr® oL —=) (19) 
rs 1 \2 1 \2 ° 
4+(aL—<) + (ol —<) 


If, starting from zero, the frequency is increased, the active component of the total impedance will in- 
crease up to a definite maximum and will then decrease somewhat, tending to the magnitude rz. The wave- 
trap suppresses all frequencies higher than the frequency at which the impedance r is the least attainable. At 
yery high frequencies, the impedance of the wave-trap decreases due to the influence of the power coil's 
natural capacitance. 


The minimum value of frequency w,, starting with which the wave-trap impedance at all frequencies will 
be greater than rg, can be determined from Eq. (19): 


ol 
OV a 


is the natural frequency of the tuned LC loop. The optimum value of the capacitance C 





(20) 





where Ww, = 


will be that for which the frequency w, is minimal. By investigating the minimum value of Formula (20), we 
get the following value for the frequency wp»: 


I. 
Woopt = “ OF Ig = WoL. (21) 
Substituting (21) in (20), we obtain 
. __ty _ L 
Wy min = “» = L orc = a . (22) 


The active component of the total impedance, as a function of frequency, has a maximum at the frequency 


If Eq. (22) is taken into account, the dependence of the wave-trap's total impedance on frequency can be 
given in the form: 


Zw 
—_—_ = 


‘ g 
a = Boe ti ecesT: (24) 





where — = = = oe is the relative frequency. 
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Fig. 7 shows the dependence on frequency of the total impedance, and its components, for a wave~trap 
based on a high-pass filter circuit. 
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5. Comparison of the Various Circuits for Tuning High-Frequency Wave-Traps 





For a comparison of the various circuits for tuning a wave-trap, Fig. 8 shows the dependence on frequency 
of the active and capacitive components of the total impedance of the wave-trap with a power coil of 0.25 milli- 
henry for various tuning designs for the frequency fp = 100 kilocycles, 
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Fig. 8. Dependence on frequency of the active and capacitive components of the total impedance 
of the wave-trap when the inductance of the power coil equals 0.25 millihenry; 1) a sharply reson- 
ant circuit; 2) the dulled circuit of Fig. 2b; 3) the broad-band tuning of the circuit of Fig. 4a with 
k = 0.6; 4) the broad-band tuning of the circuit of Fig. 4b with k, = 0.9 and ks = 0.7. The solid 
lines denote active components, the dotted lines the capacitive components. 


The curves of Fig. 8 show that, for an admissible magnitude of capacitive impedance of 250 ohms (the in- 
sertion loss with this being equal to 0.25 nepers), all the tuning designs give approximately equal trapping bands 
for the capacitive components, although at the band limits the trapping by the capacitive components of the 
two- and three-loop circuits is higher than in the circuits with one resonant loop. The width of the trapping band 
by the active components of the total impedance is significantly higher with broad-band tuning than with the 
simple resonant loops. If we set the minimum admissible wave-trap impedance at the band limits equal to 
tw min = 700 ohms and take the trapping band of the dulled circuit as being unity, then the relationships of the 
trapping bands for the different tuning schemes can be expressed in the following manner: for two-looped, 1.55; 
for three-looped, 2.0; and for a simple resonant circuit, 0.31. 


With a power coil with high inductance, when the lower boundary of the trapping band is close to the fre- 
quency w,, determined from Formula (20), the largest trapping band is given by the high-pass filter circuit, 
since in this circuit the upper limit on the band is unrestricted, 


In comparing the various wave-trap tuning designs, account should be taken, not only of the electrical 
characteristics, but also of the difficulty in physically implementing the various circuits and the question of re- 
liability in use. This latter reduces, basically, to the effectiveness of the measures for protecting the tuning ele~ 
ments from damage when there is excessive voltage in a high-voltage network. 


The magnitude of the inductance L, in the series circuit of the designs shown in Fig. 4 is inversely propor- 
tional to the magnitude of the inductance of the power coil, Therefore, when the inductance of the power coil 
is small, for example L = 0,124 millihenry, at average and at low frequencies the size of inductance L, becomes 
very large (10 millihenry and greater), Providing a large inductance with a small natural capacitance is 2 diffi- 
cult designing task; the presence of parasitic capacitance in the coil lg can materially worsen the trapping char- 
acteristics, Therefore, the use of broad-band tuning schemes for wave-traps with power coils having a 0.124 
millihenry inductance is irrational, These power coils should be used with single-frequency or two-frequency 
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tuning, allowing for the capacitive component of the total impedance. With power coil inductance of 0.25 to 

1.5 millihenry, the use of the broad-band tuning circuits of Fig. 4 is rational, Three-loop circuits give better 

characteristics but are more complicated to prepare and require finer tuning. For power coils with inductances 
greater than 1.5 millihenry, the high-pass filter circuit is generally a rational choice. 








Fig. 9. Circuits for protecting the tuning elements from overvoltage: a) for a 
two-loop wave-trap; b) for a three-loop wave-trap. 


From the point of view of protecting the tuning elements from overvoltage, the best circuit is the dulled 
circuit and its variant, the high-pass filter circuit. The explanation for this is that, if an active impedance is 
connected in series with the condenser, an excessively high-voltage wave with a sharp leading edge which may 
be applied to the wave-trap will produce a voltage increase on the condenser whose growth is retarded. With 
this, the protective discharger, parallel to the power coil, is able to continue working as the voltage on the con- 
denser builds up to a dangerous level. 


In tuning circuits in which there is no resistance in series with the condenser, in order to retard the voltage 
increase on the condenser an additional inductance coil Lg must be connected between the power coil and the 
tuning circuit, as is shown in the circuits of Fig. 9. Connecting an additional discharger P, in parallel to the 
tuning circuit provides additional protection to the circuit [4]. 
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ANNOUNCEMENT 


The National Committee of the Soviet Union for Automatic Control and the Institute of Automation and 
Remote Control are holding a conference on the theory and application of discrete automatic systems, this con- 
ference to be held September 22-26, 1958,in the Moscow Scientific Center. 


The work of the conference will fall into the following sections; 
1. Theory of pulse and digital systems. 

2. Elements and devices of pulse and digital systems. 

3, Extremal and self-adjusting systems. 

4. Optimal relay systems. 


Those wishing to participate in the work of the conference may address their enquiries for information to; 
Moscow I-53, Kalanchevskaia, d, 15a, LAT, Laboratory 7, telephone B 3-49-52, 


The Conference Organization Committee 
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